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Abstract  of  Dissertation  Presented  to  the  Graduate  Council 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 


THE  DISTRIBUTION  OF  HOTELLING'S 
GENERALIZED  T  ^ 


By 

David  Timothy  Hughes 

August,  1970 

Chairman:   Dr.  J.  G.  Saw 
Major  Department:   Statistics 

Numerous  criteria  are  available  for  testing  the 
multivariate  general  linear  hypothesis.   In  this  work  some 
properties  of  one  of  these  criteria,  T  (a  constant  times 
Hotelling's  generalized  T  ^),  are  investigated.   The  T 
statistic  is  the  sum  of  the  latent  roots  of  a  matrix  obtained 
from  random  sampling  from  a  multivariate  normal  population, 
and  its  distribution  depends  upon  the  corresponding  popu- 
lation roots  and  some  degrees  of  freedom  parameters. 

The  orthospherical  integrals  involve  expectations 
of  the  elements  of  a  matrix  having  invariant  measure  on  the 
orthogonal  group.   These  integrals  are  defined  and  a  method 
for  evaluating  them  is  presented.   We  have  derived  a  special 
representation  of  T  which,  with  the  aid  of  the  orthospherical 
integrals,  enables  us  to  present  a  method  of  evaluating  the 
moments  of  T.   The  first  two  nonnull  moments  are  obtained 
as  an  example. 


Vll 


The  exact  distribution  of  T  in  the  general  case  is 
difficult  to  deal  with  and  for  this  reason,  three  rather 
simple  approximations  to  its  distribution  are  obtained. 
All  of  these  approximations  are  obtained  by  equating  the 
moments  of  T  with  the  moments  of  the  approximating  random 
variable.   Accuracy  comparisons  of  percentage  points  in 
certain  special  cases  indicate  the  usefulness  of  tv/o  of 
these  approximations. 

A  result  by  Wald  and  Hsu  on  the  power  function  of 
the  analysis  of  variance  test  is  used  to  obtain  an  upper 
bound  on  the  integral  of  the  pov/er  function  of  T,   A 
lower  bound  on  the  integral  of  the  power  function  of  T 
is  obtained  by  using  the  properties  of  completely  monotone 
functions.   Completely  m.onotone  functions  are  also  used  to 
demonstrate  the  monotonicity  of  the  pov;er  function  of  T. 
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CHAPTER  I 

. STATEMENT  OF  THE  PROBLEM 

1.1      Introduction 

Suppose  \<e    observe  a  number  of  m-variate  column 
vectors,  ^ii-j  •   ^  ^  l,2,...,k;  j  =  l,2,...,n.  obtained  by 
random  sampling  from  each  of  k  m-variate  normal  populations 
with  common  dispersion  matrix,  V.   Assume  that  the  eler.ents 
of  the  colum.n  vectors  associated  with  population  i:   i  =  1, 
2,..,,k,  have  expectations  given  by  the  elements  of  an 
m-dim.ensional  m.ean  vector,  £.  .   It  is  desired  to  test  the 


hypothesis 


H:ui=u2=  •••  =y 


o  -"   —~  — k 

against  (1,1.1) 

H^:at  least  one  of  the  above  equalities 
does  not  hold. 

With  . 


N  =   I  n  (1.1.2) 

i  =  l  ^ 


define  x.   and  x   bv 
—1.      —..      ' 


n. 

X.  =  --   ;   X. . ,  (1.1. 5 J 

-1.   n.  .-T  -13  '   ■ 


-  1  - 


-  2  - 


and 


1   ^ 


n- 

1 


X. 


i=l  j=l 


=  1  -iJ 


(1.1.4) 


Construct  the  m  x  ni  matrix. 


S2  = 


k  ^: 

V  V 

L  L 

i = 1  j  = ; 


(^ij  -  ^i.](21ij  -  X.J',     (1.1.5) 


where  (•)'  denotes  the  transpose  of  the  vector  in  brackets 
The  matrix  S2  X':ill  be  nonsingular  v/ith  probability  one, 
whenever  X  -  k  ^  m  (see  [1]  or  [32]).   An  appropriate  test 
statistic  for  testing  (1.1.1)  is 


i=l   ^   i'     ..         1 


-  X 


) 


=  (N  -  k)  tr  SiS2-\ 


(1.1.6) 


where  Si  is  the  m  x  m  matrix  defined  by 


=  J,  "i(ii.  -  i.JS 


X   ^  ' 


(1.1.7) 


The  null  hypothesis  is  rejected  if 


\      ^    ^' 


(1.1.8) 


where  t  is  chosen  to  satisfy 


p(t 


o 


tJHj  =  a 


(1.1.9) 


3  - 


for  some  desired  level  o£  significance,  a.   The  quantity, 

T  ^,  defined  in  (1.1.6),  has  the  distribution  of  Hotelling's 

generalized  T  ^ .    ' 
^  o 

More  generally,  let  v-:i  =  l,2,*--,Vi,  be  m-dimen- 
sional,  normally  distributed  column  v^ectors  with  ni-dimen- 
sional  mean  vectors,  y_.  ,  and  common  m  x  m  dispersion 
matrix,  V.   Suppose  S2  is  an  m  x  m  central  IVishart  matrix 
(see  [56])  distributed  independently  of  vi,-««,v   with 
dispersion  matrix  V  and  degrees  of  freedom  V2  >  m  -  1. 
Then  for  testing  the  hypothesis 


H:|Ji  =  *''=ij    =0, 
o  —  -v^    — ' 

against  (1.1.10) 


H  rat  least  one  y-  f    0, 
a  — 1  '  -' 

which  is  the  canonical  form  of  the  general  linear  hypothesis, 
an  appropriate  test  statistic  is  Hotelling's  generalized 
T  ^ ,  defined  by 

V 

1 


T  2  =  V2     y    v.  'S2- V- 

1=1 


=  V2  tr  SiSz'^  ,  (1.1.11) 


where  Si  is  the  m  x  m  matrix 

V 


1 


Si  =  I      I-Ii'.  (1.1.12) 

i  =  l   ^  ^ 


It  is  noted  that  Si  defined  in  (1.1.12)  has,  in  general,  a 
noncentral  Wishart  distribution  v;ith  dispersion  matrix,  V, 
degrees  of  freedom  v^  and  m  x  :n  noncentrality  matrix 


MM '  =   y   y  •  y  •  ' 
i=l  ""-' 


(1.1.13) 


In  the  first  example  given  in  which  T^^  was  defined 
by  (1.1.6)  , 


vi  =  k  -  1 


(1.1.14) 


V2  =  N  -  k, 


(1.1.15) 


and 


k  _        _   , 

MM'  =  I      n.  (y .  -  y  ]  [y .  -  y  )  , 
i=l   ^   ^     • 


(1.1.16) 


where 


N  J,    "i^i- 
1  =  1 


(1.1.17) 


The  general  problem  just  given  is  invariant  under 
linear  transformations.  That  is,  if  instead  of  observing 
vi,--*,^'.   and  S2  ,  we  observe  l_i,''',^^      and  S  defined  by 


z .  =  Av .  :   i  =  1 , 
— 1    --1 


(1.1.18) 


and 


S  =  AS2A 


(1.1.19) 


for  any  nonsingular  n  x  n  matrix,  A,  the  problejn  remains 

the  same.   The  m-variate  column  vectors,  z  ■  :      i  •-  l,*'*,v.,, 

— i 

have  mean  vectors 

i^  =  Ay.  :   i  =  l,-",v,  (1.1.20) 

and  the  hypothesis  (1.1.10)  is  equivalent  to 


H  :Ci  =  • • •  =  ?    -  0, 
o  —'■  — v^    -' 

against  (1.1,21) 


H  : at  least  one  C •  r    0 . 
a  —1 

This  property  leads  to  consideration  of  test 
criteria  which  are  also  invariant  under  linear  transfor- 
mations.  Any  such  criterion  must  be  a  function  of  the  roots 
of  the  determinantal  equation 

|Si  -  6S2 I  =  0,  (1.1.22) 

where  Si  is  defined  in  (1.1.12).   The  distribution  of  such 
a  criterion  is'  a  function  only  of  m,  v^,  V2  and  the  non- 
zero solutions  of 

|MM'  -  XV|  =  0.  (1.1.2  5) 

(At  most,  min(m,v^)  of  Xi , • • ' ,\      are  nonzero.) 

A  number  of  test  criteria  have  been  proposed  for 
testing  (1.1.10).   These  include: 


-    6    - 


m 


1.  Likelihood   ratio:         n    (l    +    Q.]'^;  ■         (1.1.24) 

i  =  l      ^ 

2.  Roy's  largest  root:   max[8.j;  (1.1.25) 

m 

3.  Pillai's  criteria:    I      9 .  [l  +  0.)'^;      (1.1.26) 


m 
4.   Hotelling's  generalized  T  ^:v^      I    0..      (1.1.27) 

o      i=l  ^ 

In  this  paper  i\-e  investigate  some  of  the  properties  of 


T  =  T^Vv,  (1.1.28) 

for  testing  the  general  linear  hypothesis.   (It  has  become 

standard  practice  in  recent  years  to  study  the  properties 

of  T  rather  than  T  ^  itself.   We  will  continue  in  this 

o 

practice  and  steadfastly  hold  that  the  sym.bol  "T"  will 

denote  (1.1.28).   The  symbol  "T  ^"  will  be  used  onlv  when 

^       o 

specifically  discussing  the  properties  of  (1.1.11)  rather 
than  (1.1.28).   The  symbol  "T"  should  not  be  confused  with 
Hotelling's  T^  (see  equation  (1.2.1)].   Ulien  discussing 
(T^Vv^)^  we  shall  use  "(T)^"  and  reserve  "T^"  for  (1.2.1). 

1.2      Historical    Background 

In  1931,  Hotelling  [11]  suggested  a  statistic  which 
he  called  T^  as  a  generalization  of  "Student's  t"  ratio. 
Hotelling's  T^  is  in  fact  a  special  case  of  T  ^  given  in 
(1.1.11)  v;hen  v   =  1.   Using  (1.1.11),  T^  may  be  defined  as 


T'  =  ^2y'S2-'v.  (1.2.1) 

The  distribution  of  T^  was  derived  by  Hotelling  [11]  under 
the  null  hypothesis,  and,  in  1958,  Hsu  [14]  derived  the 
nonnull  distribution. 

The  statistic  T  [see  (1.1.28))  was  first  suggested 
in  1938  by  Lawley  [26]  as  a  generalization  of  Fisher's  z 
(the  familiar  F  ratio)  for  testing  significance  in  multi- 
variate analysis  of  variance.   At  that  time,  Lawley  obtained 
the  distribution  of 


°|l°g  I"'    'r..y    '^o^l-  (1-2-2' 


and  showed  that  the  limiting  distribution  of  T  ^  as  V2  tends 

•^  o 

to  infinity  is  the  Chi-Square  distribution  on  mv ^  degrees 
of  freedom.   Hsu  [15]  also  demonstrated  this  result  and 
obtained  the  null  distribution  of  T  in  integral  form.   In 
this  same  paper,  Hsu  derived  the  first  two  nonnull  moments 
of  T,  but,  as  we  shall  show,  his  second  moment  is  in  error. 

Somewhat  later,  Hotelling  [12]  suggested  T  ^  as  a 
generalization  of  his  T^  statistic  for  measuring  multi- 
variate dispersion.   In  1951,  the  null  distribution  of  T  ^ 
^  '  o 

Avas  derived  by  Hotelling  [13]  for  the  case  m  =  2  and  for 
general  v^  and  Vg  _^  2  in  terms  of  the  Gaussian  hypergeo- 
metric  function. 


Since  that  time,  numerous  persons  have  studied  the 
distribution  of  T  ^  and  T.   Pillai  [29]  obtained  the  first 
three  null  moments  of  T  and  used  these  to  obtain  an  approxi- 
mation to  its  distribution.   Using  this  approximation,  Pillai 
and  Samson  [51]  gave  upper  S'o    and  1%    points  of  T  for  m  =  2,3 
and  selected  values  of  Vi  and  va.   In  1956,  Ito  [17]  ob- 
tained an  asymptotic  formula  for  the  null  distribution  of 

T  ^  as  V2  ^  °°  and  in  1960  extended  this  result  to  cover  the 
o 

nonnull  distribution  [18].   Bagai  [2]  derived  the  null 
distribution  of  T  ^  in  terms  of  definite  integrals  for  the 
cases  m  =  2,3.   The  first  two  nonnull  moments  of  T^  were 
given  by  Mikhail  [27]  but  his  second  mioment,  although  it 
differs  from  that  found  by  Hsu,  is  also  incorrect. 

The  distributions  of  latent  roots  of  random  matrices 
v,"ere  expressed  by  James  [20],  [22]  and  Constantine  [3]  in 
terms  of  zonal  polynomials  [5],  [21].   Constantine  [4]  used 
these  results  to  obtain  the  distribution  of  T,  in  both  the 
null  and  nonnull  case,  in  series  form.   These  series  con- 
verge, however,  only  for  |T|  <  1.   In  this  same  paper, 
Constantine  derived  the  rth   moment  of  T  in  terms  of  the 
zonal  polynomials  for  Va  >  2r  +  m  -  1  (the  moments  do  not 
exist  otherwise)  and  listed  the  first  tv.'o  null  moments 
explicitly. 

In  1965,  Khatri  and  Pillai  [23]  derived  the  first 
two  moments  of  T  in  the  linear  case  (the  case  of  only  one 
nonzero  solution  of  (1.1.23))  and  the  same  authors  used 


these  results  to  approximate  the  distribution  of  T  in  the 
linear  case  [24].   They  later  generalized  this  to  consider 
any  number  of  nonzero  roots  [25].   Most  recently,  Davis  [6] 
showed  that,  in  the  null  case,  the  density  of  T  satisfies 
an  ordinary  linear  differential  equation  of  degree  m  of 
Fuchsian  type  and  obtains  an  equivalent  first  order  system. 
These  results  are  used  to  tabulate  the  upper  5%  and  1% 
points  of  T^/Vj  for  m   =    3,4  [7].   Stein  [54]  showed  that 
Hotelling's  T^  is  an  admissible  test  of 


H  :y  -  0, 

0  —   — ' 

against  (1.2.3) 

H^:y  ^  0, 

and  Ghosh  [10]  extended  this  by  showing  that  T  ^  is  admis- 
sible for  testing  the  general  linear  hypothesis.   Das  Gupta, 
Anderson  and  Mudholkar  [5]  showed  that  the  power  of  T  ^  is 
an  increasing  function  of  each  X.,  the  remaining  A.'s  being 
held  fixed,  where  X,,''',X   are  defined  in  (1.1.23). 
Mikhail  [27]  and  Srivastava  [33]  later  claimed  that  the 
power  of  T  ^  depends  upon  A  !,•••,  A   only  through 

m 

1  A-- 
i=l   " 

This  claim  is  false,  however,  since  the  second  moment  of 


T  ^  depends  upon 


m 
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Ito  [19]  gave  power  comparisons  of  T  ^  and  the  likelihood 
ratio  criteria  for  the  linear  case  and  noted  little  or  no 
difference  in  their  powers  for  large  V2.   For  the  case 
m  =  2 ,  extensive  power  comparisons  of  T  and  three  other 
test  criteria  for  the  general  linear  hypothesis  have  been 
made  by  Pillai  and  Jayachandran  [30]. 

1.2      Summary    of  Results 

Before  proceeding  with  a  discussion  of  the 
distribution  of  T,  some  preliminary  mathematical  results 
are  needed.   Therefore,  Chapter  II  is  devoted  to  the 
"orthospherical  integrals."   These  integrals  are  defined 
and  a  method  of  evaluating  them  is  presented.   The  ortho- 
spherical  integrals  up  to  order  eight  are  tabulated  at  the  • 
end  of  Chapter  II.   In  Chapter  III,  a  general  method  for 
evaluating  the  moments  of  T  is  presented  and  the  first  two 
nonnull  moments  are  found  explicitly.   Tliese  moments  are 
used  in  Chapter  IV  to  obtain  three  approximations  to  the 
distribution  of  T.   The  SI    and  11    points  of  T  are  found 
for  each  approximation,  and  accuracy  comparisons  are  made 
with  the  exact  percentage  points  and  Pillai 's  approximate 
percentage  points  for  m  =  2.   In  the  case  m  =  3,4,  compari- 
sons are  made  with  the  percentage  points  given  by  Davis  [7]. 
Upper  and  lower  bounds  are  given  for  the  integral  of  the 
power  function  of  T.   In  addition  we  offer  an  alternative 
proof  of  the  result  due  to  Das  Gupta,  Anderson  and  Mudholkar, 


11 


1.4      Notation 

Any  work  in  multivariate  analysis  is  greatly 
facilitated  by  frequent  use  of  matrix  notation,  and  it  is 
also  very  convenient  to  abbreviate  the  distributional 
properties  of  random  variables.   For  these  reasons,  certain 
notational  conventions  have  been  adopted.   Although  many  of 
these  are  standard,  they  will  be  listed  here  for  reference. 

1.  Matrices  will  be  denoted  by  capital  letters  and 
the  first  time  a  matrix  appears,  its  row  and 
column  dimensions  will  be  given.   Thus  A:(r  x  s) 
denotes  a  matrix  with  r  rov;s  and  s  columns.   The 
matrix  having  all  its  elements  zero  will  be  denoted 
by  (0).   "I"  will  be  reserved  for  the  identity 
matrix. 

2.  The  elements  of  a  matrix  will  be  denoted  by  the 
corresponding  small  letter  with  subscripts  to 
denote  their  row  and  column  position.   Thus  a. . 
denotes  the  element  in  the  ith    rov;  and  jth    column 
of  the  matrix  A.   The  symbol  (A) . .  is  equivalent 
to  a.,  and  is  sometimes  used  for  convenience. 

3.  An  underscored  small  letter  invariably  represents 
a  column   vector.   Its  row  dimension  will  be  given 
the  first  time  a  vector  appears.   Thus  x: (m)  denotes 
a  column  vector  consisting  of  m  elements.   The 
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vector  v:hich  has  all  its  elements  zero  will  be 

denoted  by  0^. 

It  is  sometimes  convenient  to  form  rov;  or  column 

vectors  from  the  rows  or  columns  of  a  matrix. 

The  sym.bol  (A)  .   will  denote  the  row    vector 

formed  from  the  ith    row  of  the  matrix  A.   The 

symbol  (A)  .  denotes  the  column   vector  formed 

from  the  ]th    column  of  the  matrix  A. 

It  is  sometimes  convenient  to  partition  a  square 

matrix  A: (r  x  r)  as 


A 


1  1 


A 


2  1 


Ai 


A 


2  2 


where  Ai i  is  the  p  x  p  matrix  formed  from  the 
first  p  rows  and  columns  A;  A12  is  the  p  x  q 
matrix  formed  from  the  first  p  rows  and  the  last 
q  columns  of  A,  p  +  q  =  r.   A21  and  A22  are  simi- 
larly defined. 

The  Kronecker  product  of  A:(r  x  r)  and  '&:{s    x    s) 
will  be  denoted  by  A  8  B  and  is  a  matrix  C:(rs  x  rs) 
with 

^(k-l)s+i, (£-l)s+j  "  ^k£^ij- 
If  X  is  a  random  variable  having  the  normal  density 

with  mean  y,  and  variance  a^,  this  will  be  denoted 


by 


-  13  - 

8.  If  X  is  a  random  variable  on  (0,=°)  with  density 

f(x)  .  {r(|)2^-^2}-'  x-lexp|4x^}. 

we  shall  abbreviate  this  by 

This  is  to  be  read  as  "x  has  the  central  Chi 
density  on  v  degrees  of  freedom." 

9.  If  X  is  a  random  variable  on  (0,=°)  with  density 

fCx)  -   I   exp(-i]f,|)^{j!r(^]2<-  =  j>':}-" 
j  =  0       "   ^         "^ 

xx^-2^^  -  lexp  (-1). 

we  shall  denote  this  by 

This  is  to  be  read  as  "x  has  the  noncentral  Chi- 
square  density  on  v  degrees  of  freedom  v;ith  non- 
centrality  parameter,  A."   If  A  =  0,  the  density 
becomes  central. 
10.   If  Xi  and  Xo  are  independently  distributed  with 


and 


'^1 


X2  -  X^  MO), 

2 


-  14  - 
then 

Vi     X2 

has  the  noncentral  F  density  on  Vi  and  V2  degi-ees 
of  freedom  and  noncentrality  parameter,  X.   This 
will  be  abbreviated 


z  -  F      (A)  . 
11.   If  Xi  and  X2  are  independently  distributed  with 

and 

X2  -  X^  '(0), 

2 

then 

X2 

u  = 


Xi   +  X2 

has  the  noncentral  Beta  density  on  Vi  and  Va 
degrees  of  freedom  with  noncentrality  parameter, 
X.      This  will  be  denoted  by 

u  -  Be(vi ,V2  ,A3  . 

12.   If  X  is  an  m-dimensional  column  vector  whose 

elements  have  a  joint  normal  density  with  mean 
vector,  jj :  (m)  and  dispersion  matrix  V:(m  >^  m)  , 
this  will  be  denoted  by 

X  -  N  (u.,V)  . 
—    m  ^'  ^ 
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13.  Let  ^i,''',^     be  independent  n-variate  colunin 
vectors  with 

X-  -  i\  fu.  ,vl  :   i  =  l,---,n. 

Then  -.--ith  X  =  (xi,---,x^_)  and  M  =  (yi,-",i^), 
the  joint  density  of  the  eloTTxents  of  X  will  be 
abbreviated 

'  X    -   N   ^    (M,V  0    1). 

14.  With  X  an  m  X  n  matrix  such  that 

X  -  N'    (M,V  8  I)  , 
mxn   '      ^  ' 

the  iTi  X  m  natrix 

S  =  XX  ' 

has  a  noncentral  I'-'ishart  distribution  with 
dispersion  natrix  V,  degrees  of  freedom  n,  and 
noncentrality  matrix  IM' .      This  will  be  abbre- 
viated 

S  ~  !^/  (V,n,MM')  . 

15.  If  S  is  an  m  X  m  symmetric  matrix  whose  m(m  +  l}/2 
mathematically  independent  elements  have  the  density 

fCSj  =  A^t^,v)]S|^    ^'       exp  [-y  tr  V   b) 

on   S  positive    definite   v;ith 

K    (V,v)    =    |2mv/2^m(m-l)/4|,^.jv/2   ;^   pfV^l-i,^-^ 
m  j=l     ~      - 
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this  will  be  abbreviated  by 

This  is  to  be  read  as  "the  ro(m  +  l)/2  mathemati- 
cally independent  elements  of  S  have  the  central 
Wishart  density  with  dispersion  matrix,  V,  and 
degrees  of  freedom  v." 

16.   The  symbols  "a  ==>  b"  and  "a  <==>  b"  are  to  be 
read  "a  implies  b"  and  "a  implies,  and  is  implied 
by,  b,"  respectively. 
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18, 


Suppose  it  is  desired  to  integrate  the  function 
f(xi,***,x  )  over  the  set 


n' 


A   =    {i'xi,''',x    1:    a.    <    x.    <   b.: 
^  ^  n-'         1   —      1   -      1 


1, • • • ,m} 


Then 


f(x)    d(x) 
A 

will   be   used   to    abbreviate 
•bi  rb 


ai 


n 


n 


f  (>^i  ,  •  •  •  ,>^j^)  dx^'  •  'dxi  . 


For  referencing  within  the  text  [•]  will  denote 
bibliographical  references,  v/hile  (•)  will  denote 
references  to  equations.   Thus  [8]  refers  to  the 
eighth  entry  in  the  bibliography,  while  (4.2.6) 
refers  to  equation  six  of  Section  2  in  Chapter  IV. 

Other  special  notational  conventions  v/ill  be 
introduced  as  they  occur  in  the  text. 


CHAPTER  II 
ORTHOSPHERICAL  INTEGRALS 

2.1       Introduction 

Let  Q,[m   x  m)  denote  the  group  of  m  x  m  orthogonal 
matrices.   Let  $  be  an  m  x  m  matrix  of  random  variables 
with  $  belonging  to  the  group  Q  (m  x  m)  ,  and  suppose  that  <l> 
has  invariant  measure  on  f2  (m.  x  m)  .   This  measure  can  be 
determined  in  the  following  manner.   Consider  any  row  of  $, 
say  the  first  rov/  for  convenience.   Then  since  $  is  orthogo- 
nal, the  elements,  <t>iy.      j  =  l,'"',m,  of  ($)  i  .  must  satisfy 

,     m 

j  =  l    -' 
=  1.  (2.1.1) 

This  is  equivalent  to  saying  that  ($)!.  determines  a  point 
on  an  m-dimensional  hypersphere  of  unit  radius.   If  $  is  to 
have  invariant  measure  on  Q(_m   x  m)  ,  then  the  measure  on 
($) 1 .  must  be 


f[($)i.]  -  s~TlT    ^^  ('J')i.(^)i.'  =  1 
m^  '^ 

=  0   otherwise, 

where 


(2.1.2) 


S^.(r)  =  2TT-^/2^^'Vr(j72)  (2.1.3) 
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is  the  surface  area  of  a  hypersphere  of  radius  r  in.  j 
dimensions . 

Now  fix  the  elements  in  the  first  row  of  i ,    and 
consider  any  other  row,  say  the  second  row  for  convenience 
Due  to  the  orthogonality  of  $,  the  elements,  'i>zA:      j  =  1, 
•••,m,  of  ('5)2.  must  satisfy 

m 

and  (2.1.4) 

m 

3  =  1    ^ 

This  implies  that  ($)2.  determines  a  point  lying  on  an 
(m  -  1) -dimensional  hypersphere  of  unit  radius  which  is 
orthogonal  to  the  vector  (=5)i..   Since  $  is  to  have  invari- 
ant measure  on  Q.[m   x  m)  ,  the  measure  on  (*)2.  conditional 
on  ($)  1  .  must  be 

f[(*)2.  I  ('t')i.]  =  s VlT    ^^  (2.1.4)  is  satisfied; 

m-1  ^  ^ 

(2.1.5) 

-    0   otherwise. 

Continuing  in  this  manner,  fix  the  elements  of 

(<f )  1  .  ,  (^)  2  .  ,  •  •  •  ,  ('5)  ■  -,       ,    and  consider  the  ith    row,  ($)  .  , 

1  - 1 ,  •  1 . 

of  f.   Since  $  is  to  be  orthogonal,  the  elements,  * . . : 

j  =  l,''',m,  of  ($) .   must  satisfy 

1  •  ' 
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m 


j=l   ^  ^^ 


=  0, 


m 


y    0, .({). .  =  0, 


(2.1.6) 


m 


3^1  *i-lJ^J 


=  0, 


and 


m 


3  =  1    -' 


This  implies  that  ($) .   determines  a  point  on  an  (m-i+1) -di- 
mensional hypersphere  of  unit  radius  which  is  orthogonal  to 

the  mutually  orthogonal  vectors  (<5)  i .,  C<^)  2  .,••',  (<^)  • 

1  -  i  ,  • 

Since  Ois  to  have  invariant  measure  on  Q.(_m   x  m)  ,  the  measure 
on  ($).  ,  conditional  on  (<*>)  1  .,  f*^)  2  .,*••,($)  •  -,       ,    is 


(*)ij  (*)l.,(*)2.,---,(*)j.j^ 


Sm-i.lfTT 


if  (2.1.6)  is  satisfied;  (2.1.7) 

=  0  otherivise. 


Hence  the  measure  on  ^    is 


m 


j=i    r  ' 


=    0   otherwise 


if  $tt>'  =  I; 


(2.1.8) 


_   -7 


0  - 


Using  (2.1.2),  we  can  v/rite  the  measure  on  $  as 


m 


£(<!>)  =   n 


r(j/2) 


if  <!>$  '  =  I; 


j  =  l   277- 

=  0   otherwise 


(2.1.9) 


In  the  derivation  of  the  moments  of  T,  v;e  shall 
need  certain  moments  and  product  moments  of  the  elements 
of  $.   To  this  end,  for  any  m  x  m  matrix,  C,  of  nonnegative 
integers,  define 


J(C)  =  E 


m   m     c.  . 

n   n  cj).  .  ^J 
i=i  j=i 


ij 


fi(mxm) 


m   m 

n   n 
i=i  j^i 


ij 


ij 


f($)  d$.    (2.1.10) 


These  are  the  "orthospherical  integrals."   We  shall  say 
that  J(C)  is  of  order  s,  where 


m   m 


i=l  ]~l      iJ 


(2.1.11) 


2.2      Some    Relationships    on    the    Orthospherical   Integrals 

Let  P  and  Q  be  fixed  m  x  m  orthogonal  matrices  and 
define  an   m  x  m  random  matrix,  0,  by 


0  =  P<|)Q, 


(2.2.1) 


21  - 


wherein  $  is  random  with  invariant  measure  on  fi (m  x  m) .   Then 
0e$^(m  X  m)  and,  since  it  is  obtained  from  $  by  an  orthogonal 
transformation,  0  also  has  invariant  measure  on  Q (m  x  m) . 
Now 


J(C)  =  E 


m   m 

n   n 
i=i  j=i 


c.  . 


IJ 


=  E 


m   m 

n   n 
i^i  j^i 


m   m 

l^    3^^  Pia*ae^3j 


c.  . 


(2.2.2) 


By  suitable  choice  of  P  and  Q,  (2.2.2)  yields  a  number  of 
interesting  relationships  on   the  orthospherical  integrals. 

In  particular,  suppose  that  P  is  obtained  from  the 
m  X  ixi  identity  matrix,  I,  by  some  permutation  of  the  rows 
or  columns  of  I,  and  that  Q  is  similarly  obtained  by  a 
possibly  different  permutation.   P  and  Q  are  therefore 
orthogonal,  and  clearly  they  have  one  and  only  one  element 
in  each  ro\v  and  column  equal  to  unity  and  all  other  elements 
zero .   That  is , 


p.   =  1 


if  a  =  i 


for  some  i 


i  ' 


=  0     if  a  ^  i-  , 


where  i.   differs  for  each  i,  and 
1  ' 


^3j  -   ^  if  B  =  j. 


for  some  j  .  , 


(2.2.3) 


=  0     if  6  M^  , 


(2.2.4) 
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where  j •   differs  for  each  j .   Now 


ill        111 

l^      I,    Pia*aB^6j 
a=l  6  =  1  ■' 


"i   'Jj 


(2.2.5) 


;o  that  (2.2.2)  can  be  written  as 


J(C)  = 


m   m 

n   n 
i=i  j=i 


i  *  i  .* 
^i  '^j 


c  . 


(2.2.6) 


with  *.  ^  .  ^  the  elements  of  0,  where  0  is  defined  m 
1  -^  J 


(2.2.1).   We  may  rewrite  (2.2.6)  as 


J(C)  =  E 


m   m 

n   n 


Xi",Jj 


IJ 


(2.2.7) 


with  c 


1  '-'j 


i=l  j=l 

^  the  elements  of  C=^  :  (m  x  m)  defined  by 

C*  =  P'CQ'.  (2.2.8) 


That  is,  C*  is  obtained  from  C  by  a  permutation  of  the  rows 
or  columns  of  C.   Since  by  definition 


J(C=^) 


m   m 

n   n  (j).  . 


1   -^3 


i=l  j=l 


ij 


(2.2.9) 


then 


J(c)  =  J[cn 


(2.2.10) 


The  permutation  matrices,  P  and  Q,  were  obtained  arbitrarily 
so  that  (2.2.10)  holds  for  all  permutations  of  C.  Thus,  for 
example,  with  m  =  4 


23  - 


J 


2    0    10' 

0    0    0    0 

=    J 

0    3    0    0 

0    0    0    0 

i 

0 

0 

0 

o" 

0 

3 

0 

0 

1 

0 

0 

2 

0 

0 

0 

0 

=  J 


3 

0 

0 

o' 

0 

0 

1 

2 

0 

0 

0 

0 

0 

0 

0 

0 

=  J 


0 

0 

0 

o" 

1 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

3 

and  so  on. 

Equation  (2.2.8)  can  be  used  to  simplify  the 
notation  by  displaying  only  the  smallest  rectangular  array 
containing  all  of  the  nonzero  elements  of  C,  along  with  an 
index,  m,  to  denote  the  dimensions  of  $  and  C.   If  C  is  the 
smallest  rectangular  array  containing  all  nonzero  elements 
of  C,  then 


J  (C)  =  J(C) 


m 


(2.2.11) 


For  example,  with  m  =  4, 


T 


'l   1   o' 

=    J 

0    0    3 

^                             ^ 

2 

0 

0 

l' 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

0 

0  o' 

0 

0 

0    3 

0 

1 

2    0 

0 

0 

0    0 

=  J 


and  so  on. 

Let  P  and  Q  be  diagonal  matrices  whose  diagonal 
elements  are  either  ±1,  so  tliat  P  and  Q  are  orthogonal. 
Then  from  (2.2.2)  and  (2.2.11), 


J    (C)    =    E 
m^  ^ 


m   m 


c.  .    c.  . 


n   n  [p.  .q. .]  -^*-  • 


m   ni 


=  J 


m 


(c)  n   n  (p.  -q. . 

i=l  j=l  ^^  JJ 


c.  . 

13 


where 


m 


c .  m 


=  J  (c)  n  p,^.  ^-n  q 


.1 


m^  'i.l  ^11 


j  =  l   JJ 


(2.2.12) 


and 


m 


1-    j£i   13 


(2.2.13) 


m 


I      c.. 
i  =  l   ^J 


Suppose  that  c.   is  an  odd  integer  for  some  ie{l,*'',m} 
For  convenience,  assume  that  Ci.  is  odd.   Setting 


Pi  1  =  -1 , 


p.  .  =  +1 
q.  .  =  +1 


i  =  2  ,  •  •  •  ,  m , 
j  =  1 ,  •  *  •  ,  m , 


(2.2.14) 


25    - 


equation    (2.2.12)    gives 


J    (C)    =    -J    CO 


(2.2.15) 


It  follows  that 


J^(^)  =  0 
m 


(2.2.16) 


whenever  any  row  or  column  of  C  sums  to  an  odd  integer, 
Thus,  for  example, 


J 


m 


and 


J 


m 


2  1 
0  3 


4  2  0 
0  0  1 


0, 


By  considering  other  special  matrices,  P  and  Q, 
many  other  relationships  can  be  obtained  on  the  orthospheri 
cal  integrals.   For  example,  choosing 


/2     /2 


(0) 


/2     /2 


and 


(0) 


(2.2.17) 
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equation  (2.2.1)  yields 


and 


1  1 


/2 


(*1  1      +      il^2  l) 


(2.2.18) 


Hence 


2  2 


/2 


((1)1  2      -      <l>22) 


m 


2    0 
0    2 


=     ff(8ii      022      ] 

=     ■j-£'((})ii      +      Cj)2l]        (012      -      (J522] 


(2.2.19) 


Expanding  (2.2.19)  and  making  use  of  (2.2.16)  gives 


2  J 


m 


'I    I' 

a  1; 


J  (2,2)  -  J 


2  0 
0  2 


(2.2.20) 


Thus  J  f^  M  can  be  evaluated  once  J  (2,2)  and  J  [^    "1  are 
known. 

Other  relationships  can  be  obtained  by  simply 
making  use  of  the  orthogonality  of  ■!> .   Some  examples  are 
presented  below. 

m 

J 


2       2 

1.  ^\i      =  (!>  1 1 


ill 

=  1   J 


m 


2  ,   2 


,      t         r'  2       ^ 

(pii  +  2.  <}'ii(t'i-i- 
j  =  2        J 


(2.2.21) 
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Taking  expectations  in  (2.2.21)  yields 


J^(2)    =    J^(4)    .    (m    -    1)    JJ2,2) 


(2.2.22) 


2.  011      (|)12         =      4>11      4:i2 


III 


I  *. 


1=1 


12 


m 


=    <t>iy    (J5i2      +      I      <fi]    4>i2    (fJ^^o    *         (2.2.23) 
1-2 


Taking  expectations  yields 


J  (2,2)  =  J  (2,4)  +  (n  -  1)  J  (2  2] 


m 


m  *■  0  2  • 


(2.2.24) 


m 


3.   0  = 


I    *li<}'2i 
j=l      J     ^ 


=  4)1  l4'2  1 


m 


j  =  l 


J  '3 


i 


m 


=  4)1 1  i^a  1   +  I     4>i  14)2 14^1  ■*2-i 
j  =  2         J   J 


(2.2.25) 


Thus 


J  (2,2)  +  (m  -  1)  J  [^  M  =  0. 


m 


(2.2.26) 


Such  relationships  greatly  simplify  the  evaluation 
of  the  Integrals,  since  they  reduce  the  number  of  integrals 
which  must  be  evaluated  directly.   Direct  evaluation  of  the 
orthospherical  integrals  can  be  tedious,  but  a  relatively 
easy  metliod  is  presented  in  the  next  section. 
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2.3      Evaluation    of   Ortho spherical    Integrals 

The  orthospherical  integrals  which  need  to  be 
evaluated  directly  may  be  evaluated  by  a  number  of  methods. 
A  polar  transformation  on  the  elements  of  ^  may  be  used 
when  the  integrals  involve  only  the  elements  in  a  single 
row  or  column  of  $.   This  method  becomes  extremely  cumber- 
some, however,  when  elem.ents  of  more  than  one  row  or  column 
are  involved.   Perhaps  the  simplest  method  is  one  v;hich 
depends  upon  the  following  representation  of  *. 

Let  X  be  an  m  X  m  matrix  of  random  variables  and 
suppose 

X  ^   \xn^^°^'^  ®  ^^  (2.3.1) 

Define  an  m  x  m  matrix  R  to  be  a  lov;er  triangular  matrix 
satisfying 

XX'  =  RR'.  (2.3.2) 

Equation  (2.3.2)  can  be  •i'.-ritten  as  a  system  of  equations 
on  the  element  of  X  and  R. 
m 

2 

1  1   J 


=  1    '^ 

m 
j  =  l 


m 


j  =  l   ^^ 
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m 


.1^  ^aj^ij  =  raxr.^  .    r.^r.^:   i  =  3,...,m, 


and 

m  m 

If  the  diagonal  elements,  r..:   i  =  l,'«',m,  are  required 
to  be  nonnegative,  then  the  equations  (2.3.3)  uniquely 
determine  the  elements  of  R.   Hence,  R  is  the  unique  lov\fer 
triangular  m  x  m  matrix  with  nonnegative  diagonal  elements 
satisfying  (2.3.2).   We  shall  call  R  the  lower  triangular 
square  root  of  XX'. 

Define  a  matrix  $ : (m  x  m)  by 

X  =  R$.  (2.3.4) 

It  is  to  be  noted  that  the  elements  of  X,  and  hence  of  R, 
are  continuous  random  variables.   Therefore,  the  diagonal 
elements,  r..:   i  =  l,''-,m,  will  be  positive  with  proba- 
bility one,  and  hence  R  will  be  nonsingular  with  probability 
one.   Then  $  is  properly  defined  by  (2.3.4)  and 

*  =  R"^X.  (2.3.5) 

Using  (2.3.2)  and  (2.3.5),  f  must  satisfy 

$$'  =  (R"^X)  (X'R''M 
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=  R"^  (XX ')R''^ 
=  R' ^ (RR'JR'"^ 
=  I,  (2.3.6) 

so  that  5  has  zero  measure  except  on  a  subset  of  il  (m  x  m)  . 

Let  f(X)  denote  the  joint  density  of  the  elements  of  X 

described  by  (2.3.1),  and  P  be  any  fixed,  orthogonal, 

(m  X  m)  matrix.   Since  f(X)  >  0  for  any  real  X,  then 

£(XP)  >  0.   From  (2.3.2),  and  the  fact  that  P  is  orthogonal, 

we  have 

RR'  =  XX' 

=  XPP'X' 

=  (XP)(XP)  '..  (2.3.7) 

Thus  the  lo-.-rer  triangular  square  root,  R,  of  XX'  is  identical 
to  that  of  XPP'X'. .  Then 

X  =  Ri  — >  XP  =  Rop.  (2.3.8) 

This  implies  that  v/-hatever  o  has  nonzero  measure,  then  so 
does  qP.   If  g(9)  is  the  measure  on    ?  then 

g(*)  >  0     for  all  GEO(m  x  m) .        (2.3.9) 

The  matrix  o  win  be  shovm  to  have  invariant 
measure  on  n (m  x   m) .   The  joint  density  of  6  and  R  will  be 
obtained,  and  it  will  be  shown  that  the  elements  of  R  are 
both  mutually  independent  and  independent  of  the  elements 
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of  *.   The  joint  density  o£  R  and  $  ■■.•.-ill  be  denoted  by 
h(R,$).   This  joint  density  can  be  obtained  using  the 
density  of  X  given  by  (2.3.1),  and  obtaining  the  Jacobian 
of  the  transformation  from  X  to  R  and  <?  defined  by  (2.3.4) 
and  (2.3.6).   Denoting  this  Jacobian  by  J(X  ^  R,9),  then 

h(R,$)  =  f(RQ)  J(X  -y   R,5),  (2.3.10) 

where  f(-)  is  the  density  of  (2.3.1).   Direct  evaluation 
of  the  Jacobian  is  cunbersor.e  and  the  method  used  here  is 
due  to  Hsu  and  reported  by  Deener  and  Olkin  [8]. 

We  seek  the  Jacobian  of  the  transformation  from  X 
to  R  and  §  defined  by 

X  =  R® 
and  (2.5.11) 

with  all  matrices  m  x  m  and  R  lov;er  triangular.   Let  oX, 
6R,  and  og  be  m  x  m  matrices  of  elem.ents,  (5X)  .  .  ,  (6R)  .  .  , 
and  C<5'5)--,  denoting  small  changes  in  the  elements  of  X, 
R,  and  ®  respectively.   Suppose  that  the  changes  5R  in  R 
and  50  in  <|)  bring  about  a  change  oX  in  X  so  that  (2.3.11) 
is  preserved.   That  is, 

X  +  5X  =  (R  +  £R) (0  +  So) 
and  (2.3.12) 

(a  +  6<5)  (<5  +  5<5)  '  =  I. 
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Expanding  equations  (2.3.12)  and  dropping  terms  of  second 
order  in  the  5's,  we  find  that 


and 


X  +  6X  =  R<5  +  R(6$)  +  (5R) 


+  ^iSQ)  '    +    (6$)$'  =  I 


(2.3.13) 


Making   use    of    (2.3.11),    we    can  v/rite    equations    (2.3.13) 


as 


6X    =    R(5$)    +    (6R)    $ 


and 


(2.3.14) 


(S^)  '    +    (6$)    $'    =    0 


Hsu  has  shown  that 


J(X  ->  R,$)  =  J(6X  -^  6R,5#)  , 


(2.3.15) 


where  J(5X  ^  6R,5<I>)  is  the  Jacobian  of  the  transformation 
from  6X  to  6R  and  5$  defined  by  equations  (2.3.14),  in 
which  R  and  <!>  are  now  considered  to  be  fixed  m  x  m.  matrices 

For  convenience,  denote  6X  by  A,  6R  by  B,  and  5$ 
by  C.   Then 


\\[ith 


and 


J(X  ^  R,$)  -  J(A  ^  B,C) , 


A  =  RC  +  Be 


C$'  +  $C'  =  0 


(2.3.16) 


(2.3.17) 


-  J  j  - 


Define  m  x  m  matrices  A*  and  C*  by 

A*  =  A$' 
and 

C*  =  C$'. 
Using  (2.5.17)  and  (2.3.18), 


(2.3.18) 


and 


A*  =  RC*  +  B 


c*  +  C*  '  =  0, 


(2.3.19) 


so  that  C"  is  an  asymetric  matrix.   Furthermore, 

J(A  -^  B,C)  =  J(A  ->  A")  J(A*  ->  B,C*)  J(C*  ^  C)  . 

(2.3.20) 

The  Jacobians,  J (A  ^  A")  and  J(C*  ^  C) ,  are  the  Jacobians 
of  the  transformations  defined  by  (2.3.18).   These  are 
orthogonal  transformations  and  hence. 


J(A  ^  A*)  =  J(C*  ->  C)  =  1 


(2.3.21) 


The  Jacobian  J  (A*  ^  B,C*)  is  the  Jacobian  of  the  transfor- 
mation defined  by  (2.3.19).   Recall  that  B  =  6R,  and  hence 
is  lower  triangular.   Direct  evaluation  of  this  Jacobian 
can  be  carried  out  rather  easily.   It  is  in  fact 

J(A*  ->  B,C") 

3[a,.=^a.,=^a..=^...a^^*a^2*. 


aTb^ii^bar^^bza' 


a„„*ai2*ai 3* • 

mm 


■im" 


•  b   *  b   • 
ml   m2 


b   *  c  1 2  *  c  1 3  *  •  •  • 
mm  '■ "       '■  "^ 


1 T,,  '^  2  3 


•m 


2m 


m- 1  ,m 


■■) 


p    «  •  •  •  (-*  A  1 

2m     m-l,m  J 


(2.3.22) 
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■where  L  is  an  m(n-l)/2  y   in  (ni  +  l)/2  matrix  which  is 
irrelevant  to  the  value  of  the  determinant  of  the  above 
matrix,  since  the  corresponding  upper  right  hand  partition 
of  the  above  r.atrix  has  all  its  elements  zero. 

The  Jacobian  J (A*  ^  B,  C-}  is  therefore 


m 


J(A*  ^  3,C")  =   u   r. 


m- 1 


i  =  l 


11 


(2.3.23) 


Equations  (2.3.16),  (2.3.20),  (2.3.21)  and  (2.3.23)  imply 


m 


J(X  -^  R,3))  =   a   r 
i  =  l 


m-i 


11 


(2.3.24) 


The  joint  density  of  R  and  <s  is  then  obtained  from  (2.3.1), 
(2.3.10)  and  (2.3.24). 

4  tr  RR'    3   r.."'-\ 
]  1  =  1   " 


2  /-. 


h(R,<5)  =  (2^)    ^-    exp|-y  tr  RK '  I  'A      r. 


on 


Sen (m  X  m) , 


(2.3.24) 


r..  >  0:   i  =  l,*'*,m, 
11  >    >  > 


and 


-°°   <    r..  <  +<»:   i  =  2,***,m;  j  =  !,•••, i-1 

This  density  factors  into 

m 


h(R,Q)  =  hi(o) 


n  h.  [r.  .1 
.  ,  1  *-  ii-* 
1=1 


m   i-1 
n    IT   h.  .  fr.  .  I 
1  =  2  j  =  l    -' 


(2.3.25) 
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m 


hi(<I>)  =   n   — iyl^     on   $£il(m  X  m)  , 

j=l  lT^'^'^ 


(2.3.26) 


h.(r..)  = 
1  ^  11-^ 


m-  i-1 


"I  -  1 


m- 1 
r. .     exp 
11      ^ 


and 


1/2 


h.  .  [x.   .1  =  (27t)'^/2 
IJ  ^  IJ  -^  1 


1  2 

—  r .  . 

2  13 


"2  ^ii' 


on  r . .  >  0  , 
11 


on   -co  <  r .  .  <  +<: 
ij 


Thus  $  has  invariant  measure  on  f2  (m  x  m)  ,  and  its  elements 
are  independent  of  the  elements  of  R.   Furthermore, 


and 


r.  .  -^  X   ■  _._T  CO)  :   i  =  1 ,  •  •  •  ,m, 
11    '^m-i  +  1  5    5  > 


r^.  -  N(0,1):   i  =  2,---,m;  j  =  l,---,i-l. 


(2.3.27) 


Using  (2.3.4)  along  M\t\v    the  distributional 
properties  of  X,  R  and  <!> ,  the  orthospherical  integrals 
can  be  evaluated.   Some  examples  follov/. 


1.   Xi  1  =  ri  i(J)i  1 
Hence , 


ffCxii^)  =  5(rii2(i)ii^) 

=  ff(rii^)5(c})ii2) 
=  E(rii2)  J  (2). 


(2.3.28) 


m 


(2.3.29) 
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Equations  (2.3.1)  and  (2.3.27)  imply  that 


and 


ff(xii^)  =  1 


^(ri  1^)  =  m. 


Substitution  into  (2.3.29)  yields 


J    (2)    =   -. 


2.      ^(xix'*) 


1  1  "^t^i  1  "*) 


=   ff(rii'^)    ff(<t>iiM 
=   ff(rn'^)    J    (4). 


m 


But 


and 


^(xii**)    =    3 


£'(rii'^)    =   in(m   +    2), 


implying  that 


(2.3.30) 


(2.3.31) 


(2.3.32) 


(2.3.33) 


J  C4)  =  ^/^  ,  TT 

m^  -^        m(m  +  2) 


3.   Xi2  =  ri  i(J) 


1  2 


(2.3.34) 
(2.3,35) 


From  (2.3.28) , 


^Cxii^xiz")  =  £'(rii'*(f)ii2(|)i2^)  ^        (2.3.36) 


and  hence 
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=  FCrii"*)  J^(2,2).        (2.3.37) 

Since  Xn    and  X12  are  independently  distributed  as  NC0,1), 

J  (2,2)  =  —.-L^     .  (2.3.38) 

The  value  of  J  (2,2)  could  have  been  obtained  from 
m 

the  relationship  given  by  (2.2.22).   We  have 

■Jni<^2)  =  J^(4)  4.  (m  -  1)  Jj^(2,2).       (2.3.39) 
Substitution  of  (2.3.31)  and  (2.3.34)  into  (2.3.39)  yields 

i=  m(m'.  2)  '    ^^    -    1)  Jm^2,2),        (2.3.40) 
so  that 

Other  orthospherical  integrals  may  be  evaluated 
in  a  similar  manner.  The  orthospherical  integrals  up  to 
order  eight,  necessary  in  order  to  obtain  the  first  four 
moments  of  T,  are  contained  in  Table  2.1. 
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TABLE  2.1 
Orthospherical  Integrals 


(The  integrals  up  to  order  eight  not  included  in  this 
table  are  equal  to  zero) 

Ai  =  m(m  +  2) 

Aa  =  mCm  +  2)  (m  +  4) 

A3  =  m(in  +  2)  (m  +  4)  (m  +  6) 

Order  2: 

J  (2)  =  i 

Order  4: 


■^m^  2^         (m-l)Ai 

J  fi  n  = I 

mW  I''      (in-l)Ai 


Order  6: 

J  (6)  =  P- 
m^  '^    A2 


J  (4,2)  =  -r- 
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TABLE    2.1 — Continued 


m'-o    a-*  (m-1)  A2 


Cm-1)A2 


J    (2,2,2)    =   ~ 


T    ra    2-)    ^       (m+1) 
mS    2>         (m-l)A2 


m 
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020 
002 


(m^+5m-2) 
(in-2}  (m-l)A2 
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m 


J 


m 


1 

1    0 
0 


]    = 


(m-l)A2 

Cin+2) 


(m-2) (m-l)A2 


(m-2) [m-TTA7 


Order    8: 


J    (8)    =    i^ 

^  A3 


J    (6,2)    =    iA 

m      '    ^         A3 


J    f  ^    ° ")    =    15  (m+5) 
m'-o    2-^  (m-l)A3 
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TABLE    2.1 — Contir.ued 


J,(4,4) 


9 

a7 


9(m+5) (m+5) 
(m-1)  (m+l)A3 


J    (4,2,2)    =    -^ 


T    f*    2-)    =    3(Tn+l) 
^m>-2    0-'  (m-DAa 


3(m+3) 
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3(m^+7m+14) 
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(m-l)A3 


9(mH-5) 
Xm-1) (m+l)A3 


J    fs    2     i>    ^  3 

mS    0    1  "^  (m-1)  A; 
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TABLE  2.1 — Continued 


5Cm+5) 

(m-1) (m+l)A3 


m 


r  3    10 

1      1      0 
0      0      2 


1    =  5(m^  +  5Tn+2) 

(ni-2)  (m-1)  (m+DAa" 


J    (2,2,2,2) 


A3 
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2       2      2       0 
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CHAPTER  III 
MOMEXTS  OF  HOTELLIXG '  S  GENER.-\LIZED  T 

2.1       Intvoduotion 


In  Chapter  I  we  defined  Hotelling's  generalized  T 


2 


(see  equation  Cl-l-H)]  as 

Vi 

T  -  =  --J,   y   V.  'Sa'H^. 
1  =  1 

=  V2  tr  SiS2-S  (3.1.1) 

where,  under  normal  assunptions  on  the  underlying  random 
variables,  Si  and  S2  are  p.  x  m  independent  random  V^ishart 
matrices  with 


S:  ~  •/^(V,vi,M:4') 
and  (3.1.2) 

Sa  -  ^'^(V,V2,(0)). 

It  is  assumed  that  V  is  positi\'e  definite.   In  view  of  the 

invariance  of  the  general  linear  hypothesis  and  T^  under 

linear  transformations  (see  Section  1.1),  we  may  simplify 

the  derivation  of  the  m.oments  of  T  ^ . 

o 

Since   V   is   positive    definite    and 

r.&r  =     y     u.u-  '  (3.1.3) 

i=l  -^-^ 
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(see  equation  (1.1.13))  is  nonnegative  definite,  there 
exists  an  m  X  m  matrix,  G,  satisfying 

GVG'  =  I 
and  (3.1.4) 

GMM'G'  =  A 

where  A  is  an  m  x  m  diagonal  matrix  v/ith  diagonal  elements 
A !,•••, A   defined  in  (1.1.23),   We  may  therefore  assume, 
without  loss  of  generalitv,  that  in  our  definition  of  T  ^ 

°  '  '  0 

given  by  (3,1.1),  the  column  vectors  v.:   i  =  l,**',\'i  and 
the  m  X  m  matrix  S2  are  independently  distributed  with 

Zi  '^  \ki'l)  '         i  =  l,-"-,Vi,        (3,1.5) 


and 


S2  -  P/^[I,V2,(0)],  (3.1.6) 


where  _£ -  :   i  -    l,'*',vi,  satisfy 


V  1 


I   C-C-  '  =  A.  (3.1.7) 

i  =  l   ^  ^ 


Write 


and 


Y  =  Wi'-'-'i-^J  r^.i.s) 


[ii,---  i^J  (3.1.9) 


so  that  Y  is  an  m  X  vi  matrix  of  random  variables  and  5  is 
an  m  X  V 1  matrix  of  means,  and 
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Y  ^  \xv/='^  ®  I)-  •    (3.1.10) 


From  (3.1.1)  we  may  write 


T  2  =  V2   y   V.  'Ss'^v 
1=1 


_  1 


=  V2  tr  Y'Sa  -Y.  (3.1.11) 

In  this  chapter,  a  method  is  developed  for  evaluating 
the  moments  of  T^,  or, more  precisely,  the  moments  of 

T  2 
T  =  -2-  (3.1.12) 

V2 

under  the  assumptions  (3.1.6)  and  (3.1.10)  with  T  '  defined 
by  (3.1.11).   The  first  two  moments,  ^[T]  and  £"[(1)^],  are 
found  explicitly.   It  is  noted  that  the  rth   moment  of  T 
exists,  provided  Va  >  2r  +  m  -  1. 

2.2      A    Method   for   Evaluating    the   Moments    of   T  ^ 


We  have 


T  2 
T  =   ° 


V2 

-  tr  Y'S2-^Y,  (3.2.1) 

with  Y  and  S2  distributed  as  in  (3.1.10)  and  (3,1.6).   Let 
di,'*',d   be  the  solutions  of  the  determinantal  equation 

IS2  -  dl|  =  0  (3.2.2) 

so  that  di,''',d   are  the  latent  roots  of  S2.   Define  D  to 
be  an  m  X  m  diagonal  matrix  with 
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D  = 


di  0   •  •  •  0 
0   d,  •  •  •  0 


0   0   •  • 


in 


(3.2.3) 


If  £.  is  an  in-dimensional  random,  normed  eigenvector  associ- 
— 1  '  •» 


ated  with  d.,  then 
1 ' 


LS2L'  =  D, 


(3.2.4) 


where  L  is  the  m  x  m  randoin  matrix  satisfvins 


L  =  (£ 


1  > 


and 


LL'  =  I, 


(3.2.5) 


(3.2.6) 


so  that  L  is  orthogonal.   L  and  D  are  uniquely  defined  by 

(3.2.4),  (3.2.5)  and  (3.2.6)  if  we  require  di  >  d2  >  ••• 

>  d_  and  the  elements  of  £^ . :   j  =  !,•••, m, of  the  first 
m  Ij    -'     '    '  ' 

row  of  L  to  be  positive.   The  joint  density  of  L  and  D 
will  be  obtained,  and  it  v/ill  be  shown  that  the  elements 
of  L  are  independent  of  the  elements  of  D.   The  Jacobian, 
J(S2  -^  L,D),  of  the  transformation  from  S2  to  L  and  D  will 
be  obtained  by  the  method  used  in  Section  2.3. 

Let  6S2,  5L,  and  6D  be  m  x  m  matrices  of  elements 


('SS2)-.,  (51)..  and  (5D)  .  .  .  representing  small  changes  in 
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the  elements  of  S2,  L,  and  D  respectively.   We  require  that 
(L  +  5L)  remain  orthogonal  and  (D  +  6D)  remain  diagonal. 
Suppose  that  the  changes  oL,  in  L,  and  oD,  in  D,  cause  a 
change  SSa  in  S2  so  that  (3.2.4)  and  (3.2.6)  are  preserved. 
Then 

S2  +  5S2  =  (L  +  5L)'(D  +  5D) (L  +  5L) 

I  =  (L  +  5L)  (L  +  5L)  '  .  (3.2.7) 

By  expanding  (3.2.7)  and  dropping  terms  of  order  greater 
than  one  in  the  5's  and  using  (3.2.4)  and  (3.2.6),  we  have 

6S2  =  L'D(5L)  +  L'(6D)L  +  (5L)'DL, 
and  (3.2.8) 

0  =  L(6L) '  +  (5L)L'. 

As  in  Section  2 . 3 

J(S2  -^  L,D)  =  J(6S2  ^  6L,  6D)        (3.2.9) 

where  J(5S2  ^  SL,  6D)  is  the  Jacobian  of  the  transformation 
from  5S2  to  5L  and  5D   defined  by  (3.2.8),  in  which  L  and  D 
are  now  assumed  to  be  fixed  m  x  m  matrices.   For  convenience, 
denote  0S2  by  A,  6L  by  B,  and  5D  by  C.   Then  from  (3.2.8) 
and  (3.2.9)  we  seek  the  Jacobian  J  (A  -^  B,C)  defined  by 

A  -  L'DB  +  L'CL  +  B'DL 
and  (3.2.10) 

0  =  LB'  +  BL'. 
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Define  m  x  m  matrices  A*  and  B*  by 

A"  =  LAL' 
and  (3.2.11) 

B*  =  BL'. 

Using  (3.2.11),  equations  (3.2.10)  may  be  written  as 

A*  =  DB'==  +  C  +  B*  'D 
and  (3.2.12) 

0   =  B*  '  +  B*, 

so  that  B=^  is  an  asymmetric  matrix.   Then 

J(A  ->  B,C)  =  J(A  -^  A*)  J(A*  ->  B*,C)  .J(B'='  -^  B)  (3.2.13) 

where  J (A  ^  A*)  and  J(B"  ^  B)  are  the  Jacobians  of  the 
transformations  defined  by  equations  (3.2.11),  and 
J(A*  ->■  B*,C)  is  the  Jacobian  of  the  transformation  defined 
by  equations  (3.2.12).   Equations  (3.2.11)  represent 
orthogonal  transformations  and  hence 

J(A  ->  A'^)  =  1 
and  (3.2.14) 

J(B*  -^  B)  =  1. 

Equation  (3.2.15)  then  reduces  to 

J(A  ^  B,C)  =  J(A'^  ^  B*  ,C)  .  (3.2.15) 
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Direct  evaluation  of  J [A-  -^   B*,C3  is  straightforward. 
Recall  that  C  =  oD  and  hence  is  diagonal,  B-  is  asymmetric 
and  A*  is  svn— etric. 


J (A*  ^  3",C)  = 


a   (.ax  l"322' 


3U 


11   -2  2 


-a^^^aai--a3i-a32^--a^.^^a^,'^ 
'^mm  b^^=^^3i=^b32---b^^-b^2'^ 


a     *) 


"'^in,m-l') 


(3.2.16) 


Laying  out  the  matrix  of  derivatives  in  tabular  form,  we 
have 
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I 

S 


T3 


CS        o        o 


o       o       o       o 


o       o 


T) 


E 


o       o       o       o 


O        13 


■^3 


ClJ         o        o 


to 

o       o       Id       o 


-o 


'O 


rt       o       o 


o       o       o       o 


nS  CD 


to 

rt        o        o 


CM 

13 
I 
to 

O  CD  O  XJ 


O  O 


rt        o       o 


o       t3       o       o 


o       o 


rt        o       o 


O  CD  O 


o       o 


■K 

eg 


O  O  CD  O 


O  O 


g3         I— I         O 

I 


o       o       o       o 


I— {  (Nl 

1—1  CSl 

o       o 


•K  4:  4: 

E           i-H           rH  CN] 

E       rvi       to  to 

u        ^  ^  ^ 


«  4: 

i-H  CM 

E        E 

/:3  -n 


I 

E 


The  Jacobian  J (A*  -  B*,C)  is  the  absolute  value  of  the 
determinant  of  the  above  matrix  and  is 


J(A*  ^  B*,c)  =  n   n  [d.  -  d.l 

l<i<j<m        ■> 


(3.2.17) 


Thus,  from  equations  (3.2.9),  (3.2.15)  and  (3.2.16), 


J(S2  ->  L,D)  =  n   n  [d.  -  d.] 

l£i<j£m    -^  J 

From  (3.1.6),  the  density  on  S2  is 


(3.2.18) 


f(S2)  = 


2^V2/2  ^  m(m-l)/4   |^'   ^j-V;  -^  1  -  j. 

j  =  l         2 


X   |S2p2-^-^/2 


exp 


■7   ^^  S2 


(3.2.19) 


on  S2  positive  definite.   The  joint  density  of  L  and  D, 
g(L,D),  is,  from  equation  (3.2.4) 


g(L,D)  =  f(L'DL)  J(S2  ->  L,D) 


(3.2.20) 


where  f(.)  is  the  density  of  (3.2.19).   Recalling  that  L 
is  orthogonal,  we  find  that 


g(L,D)  = 


2mv2/2  ^m(m-l)/4 


m 


n  r  i^'  ^  I  -  I] 


n   n  (d.-d. 


(3.2.21) 
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on   L    satisfying   LL '    -    I,    i        >    0:      j    =    l,---,m   and   dj    >    da 
>    •••    ^    "^ni   ^    ^'      "^^^    density   of    (3.2,21)    factors    into 


g(L,D)    =    gi(L)    gaCD) 


(3.2.22) 


where 


in  .         .  ,- 

gi(L)  =     n     raj/TrJ/- 


j  =  i 

if   LL'    =    I ,    £, .    >    0 
0        otherwise, 


and 


g2(D) 


^mv2/2      n      pfVa    +    1    -    1 


n     r(- 


)  r[i) 


=    1,  •  •  •  ,m 

(3.2.23) 

m/2 


n    a^^-^-"-^^^  exp(4d  )       n      n     [d.-d  ) 

j=l      J  ^      J         l<i<i<m        ^      J 


if   di    >    da    >    • • •    >    d      >    0 


=    0        other;\[ise 


(3.2.24) 


Hence  the  elements  of  L  are  independent  of  the  elements 
of  D. 

Suppose  that  K  is  any  m  x  m  random  orthogonal 
matrix  satisfying 


KSaK'  =  D, 


(3.2.25) 


Write 


K  =  [lii,k2,---,lij^] 


(3.2.26) 
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where  k.  is  an  m-dimensional ,  random,  norined  eigenvector 
associated  with  d..   Comparing  (3.2.26)  v;ith  (3.2.25),  we 
have 

Pfk.=Jl.l=i:   i  =  l,---,m 


and  (3.2.27) 

The  measure  on  K  is  obtained  from  (3.2.23)  and  (3.2.27). 
It  is  in  fact 

h(K)  -  —  n  ^[\]h'^'  if  IvK'  =  I 

2"^  3  =  1    ^ 

=  0    otherwise  (3.2.28) 


so  that  K  has  invariant  m.easure  on  Q  (m  x  m)  and  is  indepen- 
dent of  D. 

With  Y  defined  by  (3.1.10)  and  K  defined  by  (3.2.28), 
let  Z  be  an  r?.  x  Vi  matrix  of  random  variables  satisfying 

Z  =  KY.        '  (3.2.29) 

From  (3.1.10),  the  distribution  of  Z,  conditional  on  K,  is 

ZlK  -  N      (KH,I  0  I) .  (3.2.30) 

The  unconditional  density  on  Z  can  be  obtained  from  (3.2.30) 
and  is 
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f(Z)  = 


Kef2  (mxm) 

X      n   ^UX-  dK, 
j  =  l  2tt-'^ 


^  ^  .  - mv 1 / 2 
C2ttJ     '   exp 


•^  tr(Z-KH)  (Z-K5)  ' 


(3.2.31) 


It  should  be  noted  that  since  Y  and  S2  are  independent,  Y, 
K  and  D  are  mutually  independent.   This  implies  the  inde- 
pendence of  Z  and  D. 

Using  equations  (3.2.1),  (3.2.25)  and  (3.2.29),  we 
have 


T  = 


tr  Y'Sz'^Y 


=  tr(Z'K) (K'D-^K) (K'Z) 
=  tr  Z  'D- ^Z 
m   V  1  z  .  .  ^ 
1=1  i=l   1 


(3.2.32) 


Writing 


then 


U.2  =    y    z. . 
1=1     -' 


m  u. 


T  =   I 


.  ,  d. 
1  =  1   1 


(3.2.33) 


(3.2.34) 


This  representation  of  T,  along  with  the  results  of 
Chapter  II,  and  the  distributional  properties  of  Z  and  D, 
can  be  used  to  evaluate  the  moments  of  T, 
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S.S  The    First   Two   Moments    of   T 

Using  the  representation  of  T  given  in  equation 
(3.2.39)  the  first  two  moments  of  T  are 


^(T) 


m      u .  ' 

.Id- 

1  =  1      1 


m 


.r  1 


.1  'W)  n^)^ 

1  =  1  1 


and 

ff[CT)n 


m      u-  ■ 
1  =  1      1 


(3.3.1) 


m     u. 


I      -~-    II 


2  2' 

u.    u. 

1     1 


i  =  l    d 


2        r-.?     "3^ 


irj         1    J 


1  =  1  1  1?^J  -  1    j 

(3.3.2) 

From  equations  (3.2.30)  and  (3.2.33),  u.^  conditional  on 
K  is  the  sum  of  Vi  independent  normal  random  variables. 
Thus 


u .  • 
1 


K    -    X, 


(KAK') 


11 


:       i    =    1 , • • • ,m  (3.3.3) 


vjhere    from    (3.1.7) 
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For   convenience,    let 


Y-    =    (KAK').. 
1  11 

m 


j=i    ^J    ^ 


since  A  is  diagonal.   Then 


C3.3.4) 


(3.3.5) 


u. 


K  -  X^^y^]-       i  -  1,2,... 


,m. 


(3.3.6) 


Furthermore,  conditional  on  K,  u.^  and  u.^  are  independent, 

'   1       J  ^ 

Hence,  using  (3.3.6) 


H-i'] 


=  E 


K  u^2 IK 


^  (^i^) 


K        ^ 


^K"*) 


=  E 


K  U^2 


K 


E 


vi  (Vi  +  2)  +  2Yj^(vi  +  2) 


(3.3.7) 


(3.3.8) 


and 


fu.^u.n 
^    1      2    J 


E        E  [u.^u.^] 

K  u. ,u.  K   ^   J 


=  E      E       [u^]       E       (u^) 
K  u.  K   ^    u.  K   J 


E 
K 


V  1   +  V 


(Ti  ^  "^j^  ■"  ^i'^j 


(3.3.9) 
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Recalling  that  K  has  invariant  ricasure  on  Q.[m   x  m)  ,  we  may- 
use  the  results  of  Chapter  II  to  evaluate  the  moments  of 
Y-  and  the  product  moment,  E[y.y.]    where  y-  is  defined  in 


(3.3.5). 

K      ^           K 

m 
j=l      ^^      J 

=    J 

n 

m 

C2)    y    X. 
'      1=1    J 

=1   I   A., 


(3.3.10) 


E[y    ')     =    E 
K   ^      K 


m 


j^i    ^J    ^ 


m 


I    k.  .^.2  +  y  y  k./k.-^x„A. 


m 


\(4)  I      X^  .  J  (2,2)  y  Y  X^X 


and 


m(m+2) 


m 

I 
j  =  l 


X.2  + 


m (m+  2 ) 


I    I    X.X,,  (3.3.11) 


Ir^j 


Tj 


^(YiYi) 


m 


i  =  1 


m 


y    k .  ^  X 

I  _  -1    J  ^ 


m 


jl  '^i.^'^j.'^/  *   IJ,  ki/Xj.'^V. 
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m 


=  ^.t^.^'  J,  h'  ^  ^.(;  :>  ll^.h. 


m 


(3.3.12) 

Substitution  of  equations  (3.3.10),  (3.3.11)  and  (3.3.12) 
into  (3.3.7),  (3.3.8)  and  (3.3.9)  yields 


1   ^ 
*■  1  ^        m  .  "^  -,   1  ' 


(3.3.13) 


1  =  1 


£'[u.'*]  =  vi(vi  +  2)  '+  (2vi  +  4)  i  _):'   A. 


j-1 


^     m  , 

m(m+2)  .£^   i     m(m+2)   |^4   i  j 


(3.3.14) 


and 


2Vi    V    ^      2 


^  1   J  ^  m    -,   J 


m 

I      A. 


m+1 


(m-l)m(m+2J  4^4  ^i'^j 


(3.5.15) 


Upon  substitution  of  (3.3.13),  (3.3.14)  and  (3.3.15)  into 
equations  (3.3.1)  and  (3.3.2),  the  first  two  moments  of  T 
are 


Em  = 


-  j^i    ^^ 


m 
E      I       [^]  (3.3.16) 


i  =  l   "i 


and 
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E[m']  = 


vi  (vi    +    2)    + 


:vi 


m 


m 


)       X. 
1  =  1      ^ 


in 


3  V      -\    2    ^        1  V    V    -i    -v 


X    E 


m 


i=l      ^i 


m 


Vl 


y   X 


i^4   ^TdT 
irj       1    3 


(3.3.17) 


and 


Evaluation  ox  the  terns 


1=1   1 


m 


i  =  l  h' 


1,1    d-d. 
'^r:      1  1 


involve  sone  special  problems.   These  will  be  dealt  with 
separately  in  the  next  section.   'Ve  list  these  expectations 
here,  however,  in  order  to  obtain  the  moments  of  T: 


m   ^ 

.1  ^ 

1=1   1 


Vz -m-l 


(3.3.18) 


and 


m 

y  -^ 

1=1  ^1 
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y  y   1  1  ^     m(m-l) 

4  ,4  d.d.     tV2-rO  (.V2-m-l 


J 


(3.3.20) 


Substituting  these  expectations  into  equations  (3.3.16)  and 
(3.3.17),  we  obtain 


m 


ff(T)  = 


mv  1  +    y    \ . 

i  =  l   " 
v?  -HI- 1 


(3.3.21) 


and 


(v2-m)  (v2-m-l)  (v2-m-3)  £"[(1)') 


=  mvi  ^  (mv2 -in^ -2in+2)  +  2mvi(v2-l) 


m 


[2vi  (Tnv2-in2-2m-2)  +  4(v2-l))   y   X. 

i  =  l 


+  (v2-m-2) 


m 

iii  "^ 


1  2 


m 

+  2    y    x.^ 
i=i    ^ 


(3.3.22) 


Numerous  checks  are  available  on  the  inor.ents  of  T, 
The  distribution  of  T,  and  hence  its  mo.-:ients ,   is  knoxvn 
in  the  following  special  cases.   If  m  =  1, 


T 


V, 


^(X) 


X.  ^0) 


(3.3.23) 
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Putting  m  =  1  in  equations  (3.3.21)  and  (3.3.22),  the  first 
two  moments  of  T  reduce  to  the  first  t'wo  moments  of  a  ran- 
dom variable  distributed  as  in  (3.3.23).   In  the  case  of 
vi  =  1,  T  ^  reduces  to  Hotelling's  T^  defined  in  (1.2,1). 


Hence  v;hen  Vi  =  1, 


T  = 


T- 


(3.3.24) 


The  distribution  of  Hotelling's  T^  is  known  to  be  [14] 


Y   2    7TJTT 


(3.3.25) 


It  should  be  noted  that  when  Vi  =  1,  there  is  at  most  one 
nonzero  solution  of 


I A  -  XI  I  =  0 


(3.3.26) 


Then  when  Vi  =  1 


Xm^^^3 
Xv/-m.l^OT 


(3.3.27) 


Putting  Vi  =  1  in  equations  (3.3.21)  and  (3.3.22),  we  see 
that  the  first  two  moments  of  T  reduce  to  those  of  a  random 
variable  distributed  as  in  (3.3.27).   The  limiting  distri- 
bution of  T  ^,  as  V2  tends  to  infinity,  is  known  to  be 
(see  [15]  or  [26]) 


lim  T  2  - 

0      Y 
V2^°°  ITlV  1 


m 

y 

i  =  l 


X. 
1 


(3.3.28) 
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In  equation  (3.3.21),  if  we  multiply  by  Va  and  take  the 
limit  ve   obtain 


E 


lim   T 

( 


m 

mv 1  +   y   A . 
i  =  l   ^ 


(3.3.29) 


Similarly,  multiplying  (3.3.22)  by  Vo^    and  taking  the 
limit,  we  have 


lim  T  " 
o 


m 


m^vi^  +  2mvi  +  2(mvi  +  2)  I    X. 

i  =  l  ^ 


m 

i=i    ^ 


(3.3.30) 


Equations  (3.3.29)  and  (3.3.30)  are  the  first  two  moments 
of  a  random  variable  distributed  as  in  (3.3.28).   As  a 
final  check  on  the  moments  of  T,  it  is  known  that  the 
distribution  of  T  is  invariant  under  the  transformation 
(see  [1]) 


(m,Vi,V2)  ^  (vi ,m,Vi+V2-m) 


(3.3.31) 


Simple  substitution  shows  that  the  moments  of  T  given  in 
(3.3.21)  and  (3.3.22)  are  in  fact  invariant  under  this 
transformation . 

In  the  null  case  wlien  \i    =    X2 
moments  of  T  reduce  to 


^   •  • 


=  ,\   =0,  the 
m     ' 


e[t\x,   -  ...  =  X^  =  0]  = 


mv  1 


(3.3.32) 


and 
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(v2-m}  (v2-m-l)  (v2-m-3)  £"((1)21X1  = 
=  mvi^ (rnvz-m^ -2m+2)    +    2mvi(v2-l), 


IT)         J 


ra 


(3.3.33) 


2.  4      Evaluation    of   Certain    Expectations    Used 
in    Obtaining    the    Moments    of  T 


In  Section  3.3,  ive  listed  the  values  of  E 


m   , 

1  =  1   1 


and 


1,1  dTar 


in  order  to  evaluate  the  moments  of  T  (see  equations 
(3.3.18),  (3.3.19)  and  (3.3.20)].   We  derive  these  expec- 
tations here.   Recall  that  di,-**  d   are  the  latent  roots 

m 


of  S2 .   Thus 


m       , 
1  =  1   1 


(3.4.1) 


and 


m 

i^i  ^i^ 


=  tr  [5,-^]\ 


(3.4.2) 


I  I   ^^  =  (tr  Sz-']'     -    tr(S2-^)^     (3.4.3) 

IVe  shall  need  some  results  in  matrix  algebra  and 
some  well  known  properties  of  the  Wishart  random  matrix 
(see  [1]  or  [32]).   Let  A  be  an  m  x  m  nonsingular  matrix 
and  partition  A  into 
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1  1 


Ai 


A  = 


A 


2  1 


A 


(3.4.4) 


2  2 


where  A. .  is  a  p .  x 


ij 


^  X  p   matrix  with  pi  +  pa  =  m.   Then 


lA|   =  |Aii I |A22  -  A21A11-IA12 I , 
or 

=  IA22I |Aii  -  Ai2A22-'A21 I . 

Suppose  that  B  is  an  m  x  m  matrix  satisfying 


(3.4.5) 


B  =  A"^ 


(3.4.6) 


Partition  B  in  the  same  manner  as  A, 


B  = 


Bi  1   1  B 1 2 


B2 1   I  B2  2 


(3.4.7) 


Then 


and 


Bii  =  (All  -  Ai2A22-'A2i)-^ 


(3.4.8) 


that 


B22  -  (.A22  "  A21A11   A12J 


With  A  defined  and  partitioned  as  above,  suppose 


A  ^  ^/^(l,v,(0)) 


(3.4.9) 


Then  Aii  and  A22  are  independently  distributed  and 
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■"^ii  ^    ^''p-  [l.v,CO)):   i  =  1,2.        (3.4.10) 

Also,  A.,  is  distributed  independently  of  (A..    -    A.   A,   -H   1 
•^^  ^jj     ii'ii'ij-'' 

with 

^■S-j  ■  -S-i-^i"-Sj)  ^  '-''p.f^^  -  Pi,(0)).    (3.4.11) 
Furthermore,  [A^.^.  -  A^.  .  A^^- ^A.  . )  is  independent  of 


A. .A. . -^A. .  and 
ji  11    ij 


A.  .A.  .  -^A.  .  -  W      fl,p.  ,  (0)1  . 
Till   n    p  .  *■  '^1 '  *-  ■'J 


(3.4.12) 


It  is  important  to  note  that  the  V.'ishart  distribution  with 
m  =  1  on  V  degrees  of  freedom  is  identical  to  the  x^-dis- 
tribution  on  v  degrees  of  freedom.   Thus  if 


u  -  f/i(l,v,0), 


then 


(3.4.13) 


u 


X./(0) 


V\'ith  these  results,  we  can  proceed  to  evaluate 
^[tr  S2-1],  E{tr[S2-']']    and  .{  (tr  S^'^j^  -  tT[S,-')^]. 
For  convenience,  we  shall  work  with  the  corresponding 
expectations  in  terms  of  A,  where  A  is  defined  by  (3.4.9). 

d   /v{(tr  A-')'    -    tr(A-i)2} 


IVe    require    ^[tr   A'^),    E    tr(A"^j2 
Let 


an( 


B  =  A-^ 


(3.4.14) 


Then,  by  symmetry, 
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E[tT   A"^)  =  ^(tr  B) 


=  E 


m 

I   b.. 
i  =  l   ^^ 


mfi"  Cb  1  1 3 


(3.4.15) 


From  equation  (3.4.8)  with  pi  =  1  and  pa  =   m   -    1, 


bii  =  (ail  -  Ai2A22'  A21) 


-  1 A .  .  -^  -  1 


(3.4.16) 


Thus    using    (3.4.11)    and    (3.4.13) 


1  1 


'^v    -m 


7iT0T' 


and  hence 


ff(tr  A-i)    =   ^^ — ^ 
^  -^         (v-m-lj 


(3.4.17) 


(3.4.18) 


The    evaluation   of  E{(tr   X'^^    -    tr(A'^)^}    is 
similar.      By   symmetry 


^[(tr  A'M^    .    ^^.(-^-ly^^    ^   £'[(tr   B)^    -    tr(B)2) 


m 

I      b.. 
i  =  l       ^^ 


m     m 


lib.. 

L.  L  -1-1 

=11=1    ^-' 


=    E 


I    I     [b..b..     -    b..M 


2  ^. 


m(m-l)    /7(biibc2    -    bia    J. 

(3.4.19) 
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With  pi  =  2  and  pa  =  m  -  2  and  using  (3.4,8],  we  have 


biib22"  bi2   -  |Bii 


(An  -  Ai2A22'^A2i]-^  I  .     (3.4.20) 


Therefore, 

E[(tr  A-')^    -    tr(A-M')  =  ni(m-l)  E 


Al  l-Ai2A2  2-^'^2  1  I 

(3.4.21) 


Denote  (Ai i  -  Ai2A2  2'^A2i)  by 


(All  -  Ai2A22   A21)  =  C 


I 


C  1  1     C  1  2 
C2  1      C2  2 


SO  that  from  (3.4.11) 


C  -  I-/2  (I  ,v-m+2  ,  (0)) 


(3.4.22) 


(3.4.23) 


Now  from  equations  (3.4.5),  (3.4.10)  and  (3.4.11), 


C   =  cii 


c  2  2 


Cl2' 
Cl  1 


(3.4.24) 


with 
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^^    ^    X^.^.o'CO) 


and 


C2  2 


Cl2 

Cl  1 


2"\ 


^v-m+1 


(0) 


(3.4.25) 


Tliis    implies    that 


E[itv   A'M^    -    tr(A-^)2)    =   mCm-1)    E 


X 


v-m+2 


X 


m(iTi-l) 
(v-mj (v-m- 1) 


v-m+1' 


(3.4.26) 


The  evaluation  of  s[tr(A-^)^)  can  be  carried  out 
using  (3.4.26)  and  evaluating  ^(tr  A"^)^.  Again  by  sym- 
metry we  have 


^(tr  A-i)-  =  Eitr   B) 


=  E 


=    E 


m 

I   b.. 
i  =  l   ^^ 


m 


I   b.,2  +  I    I    b-.b.. 
i=l   ^^     i^j   ^^  JJ 


=  m£'(bii2)  +  m(m-l)  F(biib22).     (3.4.27) 


From  (3.4.17), 


^(bir')  = 


(v-m-lj (v-m-5) 


(5.4.28) 


To  evaluate  £'(biib22)  consider  (3.4.7)  and  (3.4.8)  with 
pi  =  2  and  p2  =  m  -  2.   We  have 
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1  1  - 


bi  1   b  1 2 

b  2  1    b  2  2 


=   (All   -  Ai2A22-^'\2l)"^ 


=  c 


-  1 


(3.4.29) 


where  C  is  defined  in  (5,4.22)  and  distributed  as  in 
(3.4.23).   Another  application  of  (3.4.8)  yields 


and 


11- 


Cl  1 


Cl2 
C2  2 


_  1 


(3.4.30) 


Thus 


'2  2   - 


biib 


2  2 


C22 


C12 


1  1 


-  1 


C 1  1 C2  2 


C12 


C 1  1 C2  2 


-  1 


(3.4.31) 


The  remarks  preceding  (3.4.10)  along  with  (3.4.23)  imply 
the  independence  of  Cii  and  C22.   Also,  the  remarks  pre- 


ceding (3.4.11)  imply  the  independence  of  Cii  and  c 


2  2 


C12 


Ci  1 


C12 


C 1  1 C2  2 


•)  .   Similarly,  C2  2 


Hence  Cii  is  independent  of  [l  - 

2 

is  independent  of  [l  -  -^^~? ]    and  hence,  so  is  the  product 

Ci  1C2  2-^  '  ^ 

C11C22.      This    im.plies    that 
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^'(bi  ib2  2)  =  E   - — 


f   1   ] 

\    "•     ' 

F. 

E 

UiiJ 

[C22J 

Cl  2 


2   1 


-  2 


C 1  1 C2  2 


[  ^  1 

E 

1 

E 

[ci  1  J 

[C22J 

Cl  1 


>,  2 


Cl  1 


Cl  2 


C2  2 


(3.4.32) 


The  expectations  E 


r     ^      \ 


ICiiJ 


and 


.C22J 


can  be  evaluated  using 


(3.4.10),  (3.4.13)  and  (3.4.23),  yielding 


f       ^        \ 


and 


CiiJ 


f  1  1 


C2  2 


1 

v-m 


v-m' 


(3.4.33) 


To  evaluate 


Cl  1 


Cl  1 


C12 


C2  2 


we  need  to  make  use  of  (3.4.12).   Let 


Ui  =  Cl  1 


Cl  2 

C2  2 


and 


U2   = 


C12 
C2  2 


(3.4.34) 


From  (3.4.12)  and  (3.4.23),  Uj  and  U2  are  independently 
distributed  as 


and 


Therefore 
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^^  ^  X^.n^./CO) 


U2  -  Xi'(O) 


Ci  1 


Cl 1   -   Cl 2 


Ui   +  U2 


2  2 


(3.4.35) 


BeCl,v-m+lT  • 


(3.4.36) 


Equation  (3.4.36)  yields 


Cl  1 


Cl  1 


Cl  2 


C2  2 


(v-m) (v-m- 2) 
(v-m-1) (v-m-3) ' 


(3.4.37) 


Substitution   of    (3.4.33)    and    (3.4.37)    into    (3.4.32)    gives 


T?  r\.      u      -^    _  (v-m-2) 

i?lbiib22)    -   -[^-^rj-^v-m-DTAT^iToT' 


(3.4.38) 


and  we    obtain    from    (3.4.27),    (3.4.28)    and    (3.4.38) 


E{tr  A'') 


1^2      _ 


in  (mv-m^ -2m+2) 
'(v-m)  (v-m-1)  (v-m-3) 


(3.4.39) 


Finally  substituting  (3.4.39)  into  (3.4.26)  and  solving  for 
£'(tr(A' ^)  ^),  we  find  that 


n  (v  - 1 ) 


ntri^-^y]    =   ■(,-„)  ,v-m-ij  (v-m-3)-      (2. 4. 40) 
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Putting  A  =  S2  and  replacing  v  by  V2  in  equations  (3.4.18), 
(3.4.26)  and  (3.4.40)  we  obtain  the  equations 


^(''^  S,-')  =  ^^  (3.4.41) 


and 


F((tr  S.-M^  -  tr(S,-")^)  -   T,,':LTill-m-n        f^-"-"' 


Xv2  -ni)  (.V2  -ra-1  J 
which  are  the  same  as  equations  (3.3.18),  (3.3.19)  and 
(3.3.20) . 

3. 5      Summary 

The  representation  of  T  given  in  (3.2.34), 

m   u-  ^ 
1  =  1   1 

can  be  used  to  evaluate  higher  order  moments  of  T.   The 
expectations  involving  the  elements  Ui^,*'',u^,  can  be 
evaluated  using  (3.3.6)  and  the  orthospherical  integrals. 
The  orthospherical  integrals  necessary  for  the  first  four 
moments  of  T  are  tabulated  at  the  end  of  Chapter  II.   The 
expectations  involving  elements  of  D  can  be  found  using 
simple  properties  of  the  Wishart  distribution  similar  to 
those  in  Section  3.4. 
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It  should  be  noted  that  the  second  moment  of  T 
depends  upon 

m 

Ji  '^'- 

This  refutes  the  claim  by  Mikhail  [27]  and  Srivastava 
[33]  that  the  power  of  T  depends  upon  \i,---,\      only 


through 


m 
i  =  l   ^ 


CHAPTER  IV 

APPROXIMATIONS  TO  THE  DISTRIBUTION  OF 

HOTELLING'S  GENERALIZED  T  ^ 

o 

4.1      Introduction 

It  Avas  pointed  out  in  Chapter  III  that  the 
distribution  of  Hotelling's  generalized  T^  is  well  known 
in  certain  special  instances.   These  are  the  cases  when 
m  =  1  for  all  vi  and  V2  ;  when  vi  =  1  for  all  m  and  V2  _^  m; 
and  when  V2  goes  to  infinity  for  all  m  and  vi.   In  addition 
to  these  special  cases,  Hotelling  [13]  derived  the  null 
distribution  of  T^^  for  m  =  2,  all  vi  and  V2  ^  2.   With 
the  exceptions  of  these  few  instances,  the  distribution 
of  Hotelling's  generalized  T^  is  customarily  in  series 
form  in  terms  of  the  zonal  polynomials  as  described  by 
James  [21]  and  Constantino  [3].   Constantine  [4]  gives  a 
series  representation  of  the  distribution  of  T.   His  series 
is  convergent,  however,  only  when  |T|  <  1.   In  1968,  Davis 
[6]  showed  that  the  null  distribution  of  T  satisfies  a 
linear  homogeneous  differential  equation  of  order  m,  and 
that  Constantino's  series  reduces  to  the  solution  of  this 
equation  in  the  case  JTJ  <  1.   Davis  [7]  uses  this  result 
to  obtain  the  upper  5?,  and  1%   points  of  T  -^/vi  for  m  =  3,4 
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and  various  values  o£  vi  and  va  v.'hich  he  claims  are 
generally  accurate  to  five  significant  figures. 

In  1954,  Pillai  [29]  suggested  that  the  distribution 
of  T  be  approximated  by 

^  m(.V2-m-lj+2   mvj  ,m(v2-m-l)+2'  (.4.1.13 

It  should  be  pointed  out  that  this  approximation  becomes 
exact  for  m  =  1  as  well  as  for  the  limiting  case  when  V2 
tends  to  infinity.   However,  the  approximation  is  not 
exact  in  the  case  Vi  ^  1,  nor  is  it  invariant  under  the 
transformation 

Cra,Vi,V2)  ^  (vi  ,m,Vi+V2-m)  . 

In  this  chapter  three  approximations  to  the 
distribution  of  T  are  considered.   Upper  5%  and  1%  points 
are  found  for  each  of  the  approximations  in  the  cases 
m  =  2,3,  and  4  and  some  accuracy  comparisons  are  made. 

4.2      Approximations    to    the    Distribution    and 
Comparison    of  Tercentage    Points    of   T 

In  view  of  the  known  distribution  of  T  ^  for  the 

o 

special  cases  mentioned  in  Section  4.1,  a  number  of  approxi- 
mations to  its  distribution  are  proposed  for  the  general 
situation.   Of  course  we  would  like  the  approximation  to 
be  a  simple  one  in  terms  of  tabulated  functions  and  to 
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satisfy  as  many  known  results  as  possible.   Naturally,  we 
would  want  any  approximation  to  be  sufficiently  accurate 
for  practical  purposes.   For  these  reasons,  the  following 
three  approximations  to  the  distribution  of  T  are  considered 

Approximation    1 


T  •  ax^,  (y) 


(4.2.1) 


for  some  a,  v,  and  y^ 


Approximation    2 


T-^aF         ,t(y) 
mv  1  ,V2  -m  +  1  ^  '  '^ 


for  some  a  and  y;  and 


Approximation    2 


(4.2.2.) 


T  *  aF 


V ,V2 -m+1 


.l(T) 


(4.2.3) 


for  some  a,  v,  and  y. 

Approximation  1  was  chosen  in  view  of  the  fact 
that  the  limiting  distribution  of  T^  is 


^mv  ■ 


m 


.k'^ 


I 


It  was  suspected  that  this  approxim.ation  would  be  particu- 
larly suited  for  large  Va-   Approximation  2  was  chosen 
since  in  the  cases  m  =  1  or  vi  =1,  the  distribution  of  T 
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is  given  by  (3.3.23)  and  (5.3.27),  respectively.   It  was 
thought  that  approximation  2  would  be  fairly  accurate 
overall.   The  third  approximation  was  chosen  for  reasons 
similar  to,  and  in  hopes  of  improving,  approximation  2. 

In  approximations  1  and  3,  the  first  two  null 
moments  of  the  approximating  distribution  were  equated  to 
the  first  two  null  moments  of  T  given  in  (5.3.32)  and 
(3.3.33)  in  order  to  determine  the  unknown  multipliers, 
a,  and  unknown  degrees  of  freedom,  v.   In  approximation  2, 
the  first  moment   of  the  approximating  distribution  was 
equated  to  the  first  null  moment  of  T  given  in  (3.3.32)  in 
order  to  determine  the  unknown  multiplier,  a.   The  first 
nonnull  moment  of  each  of  the  approximating  distributions 
was  then  equated  with  the  first  nonnull  moment  of  T  given 
in  (3.3.21)  to  obtain  the  unknown  noncentrality  parameters, 
y.   Derivation  of  these  unknown   multipliers,  degrees  of 
freedom,  and  noncentrality  parameters  follow. 
Let  u  be  a  random  variable  such  that 


u  .  ax^^(y).  ^4.2.4) 

The  first  two  null  moments  of  u  are  known  to  be 

ff(u|Y=0)  =  av  (4.2.5) 

and 

E(unY-O)  =  a^'(v  +  2).  (4.2.6) 
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Equating  these  two  expressions  with  (3.3.52)  and  (3.3.33) 
we  obtain 


and 


a  =   -   (v2-l)(vi^V2-m-l) 

(v2-m) (v2-m-l) (v2-m-3)  l^././j 


_  mvi  (vg-m)  (v2-in-3)  ,.  ^  „. 

"^  "  I{;7^T)7vi+V2-m-l)-  (4.2.8) 


The  first  nonnull  moment  of  u  is  given  by 

E(u\y)    =    a(v  +  y) .  (4.2.9) 

Equating  this  with  (3.3.21)  and  using  (4.2.7)  and  (4.2.8) 
yields 


m 

I 
i  =  l 

m 


_   (v2-m) (v2-m-3)     r   , 
^  "  (V2-1)  (vi+V2-m-l)  ,^,  '^i 


^    ^   ^  (4.2.10) 


mv  1  .  ■^  T   1 ' 
1  =  1 


where  the  value  of  v  is  given  by  (4.2.8).   Thus  approxi- 
mation 1  becomes 

T  ;  ax./(Y),  (4.2.11) 

where  a,  v,  and  y  a^e  given  by  (4,2.7),  (4.2.8)  and 
(4.2.10)  respectively. 

If  u   is  a  random  variable  such  that 


"  ~  «''«,, V. -..!«•  f^-^-") 
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then  the  first  null  moment  of  u  is 


£'Cu|y  =  0)  =  a 


(v2 -m+1) 
Iv'2-m-lJ 


(4.2.13) 


Tlie  evaluation  of  a  by  equating  (4.2.13)  with  (3.3.32) 
yields 


a 


mv  1 
V2 -m+I' 


(4.2.14) 


Using  (3.3.21),  (4.2.14)  and  the  first  nonnull  moment  of  u, 


„ --  I  -.      (mvi  +y)  (v2  -ra+1) 

£:(u  y)  =  a  -^ r  .  \^  1  ^ ~ 

^  '  '  '^       mvi  (V2-m-l) 


(4.2.15) 


we  obtain 


m 


i=i  ^ 


Approximation  2  then  becomes 


T    • 


mvi 


.    V2  ~iTi+l      mVjjVj-m-l 


m 
i  =  l   ^ 


(4.2.16) 


(4.2.17) 


If  u   is  a  random  variable  with 


u  ^  aF     „_^T  (y)  , 
V  , v->  -m+1  ^  '  -'  ' 


(4.2.18) 


the  first  two  null  moments  and  the  first  nonnull  moment  of 
u  are  respectively 


F(uly  =  0)  =  a  -P^^^^-^, 
^  '  '  '  (v  2  -  m  - 1 )  ' 


(4.2.19) 


and 
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Eiu^  |y=0) 


a 


2   (v  +  2)  (vz-m^-l)^ 


(4.1.20) 


E[u\y)    =    a 


(v+y)  Cv2-m+l) 
V ( V 2  - m-1) 


(4.2.21) 


A  procedure  similar  to  that  used  for  approximation  1  yields 
the  following  values  for  a,  v,  and  y- 


a  = 


m.v  1 
V2-ni  +  l ' 


(4.2.22) 


and 


V  = 


mvi (v2 -m) 
V  1  +V2 -mv 1  - 1 ' 


V2  -m 


m 


Vi+V2-mvi-l  -t.-,      1 


V 


m 


mvi 


I   A.. 


i  =  l 


(4.2.23) 


(4.2.24) 


It  should  be  pointed  out  that  the  value  of  v  given  in 
(4.2.23)  may  be  negative  when  m  and  Vj  are  large,  relative 
to  V2  .   We  Avill  agree  to  use  approximation  3,  therefore, 
only  when  the  value  of  v  is  positive.   That  is,  v/hen 
V2  >  (n^-l)  Vi+1.   Thus,  approximation  3  becomes 


mvi    p 
V2  -m+T  ■"  V  ,V2  -m+1 


mv  i 


m 

I 
i  =  l 


(4.2.25) 
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provided   v  >  0,   where 

V  =  r-^-^ — ~ ^.  (4.2.26) 

V 1+V2  -mvi  -1  ^       ■' 

Approximations  2  and  3  become  exact  in  the  cases 
m  =  1  for  all  Vi  and  V2 ;  Vi  =1  for  all  m  and  V2 ;  and  in 
the  limit  as  V2  goes  to  infinity.   They  are  also  invariant 
under  the  transformation 

(m,Vi,V2)  -^    (vi  ,m,Vi+V2-Tn)  . 

Approximation  1  is  also  invariant  under  this  transformation 
and  becomes  exact  in  the  limit.  It  is  obviously  not  exact, 
however,  when  m=  1  orvi  =  1. 

The  following  tables  give  the  upper  5%  and  II 
points  of  T  obtained  from  each  of  the  three  approximations 
considered  in  this  chapter,  as  v/ell  as  some  accuracy  com- 
parisons.  In  the  case  m  =  2,  our  approximate  percentage 
points  are  compared  with  the  exact  percentage  points  ob- 
tained from  the  distribution  derived  by  Hotelling  [13]  as 
well  as  those  found  using  Pillai's  approximation  [29]. 
For  m  =  3  and  4,  comparisons  are  made  with  Pillai's 
approximate  percentage  points  and  those  given  by  Davis  [7]. 
Hotelling's  exact  percentage  points  are  reprinted  from 
[30] .   Percentage  points  for  approximation  1  were  obtained 
using  Percentage  Points  of  the  x^ "distribution  [28].   Th( 
percentage  points  for  approximations  2  and  3,  as  well  as 


le 


S3 


TABLE    4.1 


Upper  5 

%    points 

o£  T 

(m  =  2) 

V2 

Vi 

13 

23 

33 

43 

53 

63 

83 

125 

E 

1.4508 

.6807 

.4429 

.3280 

.2603 

.2158 

.1607 

.1064 

P 

1.3966 

.6683 

.4385 

.3261 

.2596 

.2156 

.1610 

.1069 

3   1 

1.4994 

.6861 

.4446 

.3286 

.2606 

.2159 

.1607 

.1064 

2 

1.4979 

.6952 

.4524 

.3341 

.2548 

.2192 

.1631 

.1078 

3 

1.4473 

.6820 

.4447 

.5297 

.2619 

.2172 

.1619 

.1073 

£ 

2.1771 

1 

.0071 

.6517 

.4812 

.3813 

,3157 

.2348 

.1552 

P 

2.0893 

9849 

.6424 

.4763 

.3784 

.3138 

.2340 

.1550 

5   1 

2.2284 

1 

0097 

.6515 

.4807 

.3808 

.3153 

.2345 

.1551 

2 

2.2945 

1 

0439 

.6703 

.4924 

.3889 

.3212 

.2381 

.1569 

3 

2.1741 

1 

0073 

.6525 

.4821 

.3821 

.3164 

.2354 

.1557 

E 

2.8802 

1 

3197 

.8508 

.6269 

.4960 

.4103 

.3049 

.2013 

P 

2.7658 

1 

2907 

.8384 

.6202 

.4920 

.4076 

.3035 

.2008 

7   1 

2.9324 

1 

3193 

.8488 

.6253 

.4949 

.4094 

.3043 

.2010 

2 

3.0819 

1 

3829 

.8817 

.6454 

.5084 

.4192 

.3101 

.2058 

3 

2.8828 

1 

3212 

.8521 

.6279 

.4970 

.4111 

.3055 

.2017 

E 

4.9409 

2 

2260 

1 

.4243 

1.0451 

.8247 

.6809 

.5046 

.3323 

P 

4.7566 

2 

1806 

1 

.4050 

1.0346 

.8182 

.6765 

.5023 

.3314 

13   1 

4.9920 

2 

2157 

1 

.4169 

1.0402 

.8213 

.6784 

.5032 

.3317 

2 

5.4135 

2 

3783 

1 

.4998 

1.0903 

.8549 

.7024 

.5173 

.3383 

3 

•    •    • 

2 

2342 

1 

.4288 

1.0480 

.8268 

.6825 

.5056 

.3329 

-  84  - 


TABLE  4.2 

Uoper  II  points  of  T 
Cm  =  2) 


V2 


Vl 

13 

23 

33 

43 

55 

65 

85 

123 

E 

2.2704 

.9850 

.6259 

.4558 

.5589 

.2959 

.2189 

.1440 

P 

2.0541 

.9350 

.6029 

.4445 

.5519 

.2912 

.2165 

.1430 

3   1 

2.1429 

.9740 

.6068 

.4460 

.5525 

.2914 

.2164 

.1428 

2 

2.4103 

1.0250 

.6458 

.4675 

.5666 

.2820 

.2222 

.1456 

3 

2.50  2  5 

.9953 

.6295 

.4593 

.5613 

.2977 

.2201 

.1447 

E 

3.2642 

1.3895 

.8759 

.6361 

.4998 

.4114 

.5059 

.1995 

P 

2.9624 

1.5218 

.8456 

.6208 

.4902 

.4050 

.5004 

.1980 

5   1 

3.0447 

1.5287 

.8467 

.6209 

.4900 

.4947 

.5001 

.1978 

9 

3.5800 

1.4807 

.9167 

.6609 

.5159 

.4228 

.5104 

.2025 

3 

3.3535 

1.4100 

.8840 

.6422 

.5038 

.4145 

.3056 

.2005 

E 

4.2220 

1.7736 

1 

.1099 

.8058 

.6320 

.5197 

.5855 

.2515 

P 

3.8445 

1.6922 

1 

.0764 

.7878 

.6208 

.5121 

.5792 

.2495 

7   1 

3.9124 

1.6911 

1 

.0755 

.7855 

.6191 

.5108 

.5783 

.2491 

9 

^ 

4.7468 

1.9223 

1 

.1796 

.8461 

.6583 

.5582 

.5959 

.2561 

3 

4.5327 

1.8065 

1 

.1258 

.8152 

.6582 

.5241 

.5859 

.2525 

E 

7.0198 

2.8804 

1 

.7842 

1 

.2879 

1.0065 

.8255 

.6069 

.5964 

P 

6.4298 

2.7618 

1 

.7373 

1 

.2634 

.9915 

.8155 

.6015 

.5942 

13   1 

6.4443 

2.7359 

1 

.7208 

1 

.2530 

.9844 

.8104 

.5985 

.5928 

2 

8.1696 

3.2149 

1 

.9432 

1 

.3809 

1.0675 

.8687 

.6517 

.4078 

3 

•    •    • 

2.9509 

1 

.8182 

1 

.3079 

1.0196 

.8549 

.6125 

.5989 

-  8  5 


TABLE  4.3 

Upper  5%  points  of  T 
Cm  =  3) 


V2 

Vl 

30 

40 

50 

60 

80 

100 

200 

D 

1.0385 

.7348 

.5682 

.4631 

.3380 

.2661 

.1289 

P 

1.0097 

.7201 

.5595 

.4574 

.3351 

.2644 

.1286 

5 

1 

1.0361 

.7329 

.5668 

.4621 

.3374 

.2657 

.1288 

1.0936 

.7655 

.5897 

.4768 

.3458 

.2712 

.1303 

3 

1.0407 

.7364 

.5695 

.4642 

.5388 

.2667 

.1292 

D 

1.8477 

1.3020 

1.0045 

.8174 

.5954 

.4682 

.2263 

p 

1.8022 

1.2791 

.9908 

.8084 

.5907 

.4653 

.2257 

10 

1 

1.8379 

1.2960 

1.0005 

.8146 

.5938 

.4672 

.2260 

2 

2.0026 

1.3892 

1.0605 

.8564 

.6175 

.4824 

.2299 

3 

1.8579 

1.3081 

1.0085 

.8203 

.5972 

.4694 

.2266 

D 

2.6327 

1.8503 

1.4251 

1 

.1584 

.8425 

.6619 

.3193 

P 

2.5731 

1.8206 

1.4075 

1 

.1468 

.8365 

.6583 

.3185 

15 

1 

2.6157 

1.8402 

1.4181 

1 

.1531 

.8394 

.6599 

.3189 

2 

2.8997 

2.0022 

1.5234 

1 

.2273 

.8819 

.6874 

.3259 

3 

.  .  . 

1.8616 

1.4325 

1 

.1636 

.8456 

.6639 

.3199 
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TABLE  4.4 

Upper  1%  points  of  T 
Cm  =  3) 


V2 

Vi 

30 

40 

50 

60 

80 

100 

200 

D 

1.3504 

.9589 

.7189 

.5824 

.4219 

.3305 

.1589 

P 

1.2773 

.9024 

.6973 

.5680 

.4143 

.3261 

.1578 

5    1 

1.3041 

.9152 

.7047 

.5729 

.4168 

.3275 

.1582 

2 

1.4749 

1.0061 

.7611 

.6112 

.4379 

.3408 

.1614 

3 

1.3785 

.9544 

.7288 

.5892 

.4257 

.3332 

.1595 

D 

2.2913 

1.5848 

1.2101 

.9783 

.7070 

.5534 

.2652 

P 

2.1863 

1.5341 

1.1805 

.9590 

.6970 

.5474 

.2638 

10    1 

2.1853 

1.5442 

1.1859 

.9623 

.6986 

.5482 

.2640 

2 

2.6151 

1.7620 

1.3219 

1 

.0053 

.7499 

.5808 

.2719 

3 

2.3590 

1.6221 

1.2337 

.9945 

.7161 

.5592 

.2666 

D 

3.1977 

2.2039 

1.6791 

1 

.3555 

.9778 

.7645 

.3654 

P 

3. 0638 

2.1406 

1.6427 

1 

.3320 

.9658 

.7572 

.3638 

15    1 

3.0826 

2.1469 

1.6447 

1 

.3322 

.9653 

.7568 

.3637 

2 

3.7372 

2.5019 

1.8684 

1 

.4865 

1.0514 

.8116 

.3773 

3 

•    •    • 

2.2629 

1.7163 

1 

.3811 

.9920 

.7736 

.3677 
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TABLE  4.5 

Upper  5"5  points  o£  T 
(m  =  4) 


V2 


Vi 

30 

40 

50 

60 

80 

100 

200 

D 

1.3638 

.9524 

.7313 

.5934 

.4308 

.3381 

.1629 

P 

1.3140 

.9274 

.7165 

.5836 

.4257 

.3350 

.1622 

5 

1 

1.3589 

.9491 

.7290 

.5918 

.4299 

.3375 

.1627 

2 

1.4621 

1.0062 

.7653 

.6168 

.4439 

.3465 

.1650 

3 

1.3689 

.9556 

.7335 

.5951 

.4319 

.3389 

.1631 

D 

2.4501 

1.7053 

1.3070 

1.0593 

.7679 

.6022 

.2896 

P 

2.3732 

1.6676 

1.2847 

1.0447 

.7603 

.5975 

.2886 

10 

1 

2.4357 

1.6971 

1.3017 

1.0  5  53 

.7654 

.6005 

.2891 

2 

2.7213 

1.8562 

1.4033 

1.1261 

.8056 

.6264 

.2957 

3 

.  .  . 

1.7166 

1.3143 

1.0644 

.7709 

.6041 

.2901 

D 

3.5074 

2.4358 

1.8643 

1.5095 

1 

.0930 

.8565 

.4113 

P 

3.4069 

2.3871 

1.8359 

1.4910 

1 

.0834 

.8507 

.4100 

15 

1 

3.4830 

2.4212 

1.8560 

1.5034 

1 

.0896 

.8545 

.4107 

2 

3.9696 

2.6961 

2.0518 

1.6266 

1 

.1597 

.8996 

.4223 

3 

.  .  . 

.  .  . 

1.8773 

1.5186 

1 

.0982 

.8599 

.4121 
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TABLE  4.6 

Upper  II  points  of  T 
(m  =  4) 


V2 


Vi 

30 

40 

50 

60 

80 

100 

200 

D 

1.7348 

1.1892 

.9040 

.7289 

.5253 

.4106 

.1962 

P 

1.6178 

1.1324 

.8706 

.7071 

.5139 

.4035 

.1946 

5 

1 

1.6745 

1.1595 

.8865 

.7174 

.5193 

.4069 

.1953 

2 

1.9498 

1.5036 

.9749 

.7773 

.5519 

.4273 

.2003 

3 

1.7819 

1.2147 

.9200 

.7399 

.5314 

.4144 

.1972 

D 

2.9810 

2.0356 

1.5441 

1.2434 

.8946 

.6985 

.3332 

P 

2.8126 

1.9564 

1.4985 

1.2139 

.8795 

.6894 

.3312 

10 

1 

2.8785 

1.9865 

1.5155 

1.2244 

.8843 

.6920 

.5316 

2 

3.5341 

2.3341 

1.7301 

1.3716 

.9657 

.7437 

.3444 

3 

•    •    • 

2.0957 

1.5816 

1.2690 

.9087 

.7075 

.3354 

D 

4.1858 

2.8501 

2.1581 

1.7359 

1.2471 

.9729 

.4632 

P 

3.9706 

2.7510 

2.1019 

1.6999 

1.2289 

.9619 

.4608 

15 

1 

4.0413 

2.7801 

2.1188 

1.7097 

1.2335 

.9648 

.4613 

2 

5.1018 

3.3493 

2.4732 

1.9531 

1.3686 

1.0506 

.4827 

3 

.  .  . 

•    •    • 

2.2170 

1.7760 

1.2692 

.9869 

.4666 
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Pillai's  approxiination,  ivere  obtained  from  Tables  of  the 
Incomplete  Beta  Distribution  [28].   There  are  five  entries 
in  each  cell.   For  m  =  2,  these  correspond  to  the  exact 
percentage  points  and  the  approximate  percentage  points 
found  from  Pillai's  approximation  and  approxim.ations  1, 
2  and  3  of  this  chapter.   These  are  designated  by  E,  P,  1, 
2,  and  3  respectively.   For  m  =  3  and  4  the  exact  per- 
centage points  are  replaced  by  those  given  by  Davis  which 
are  designated  by  D.   Davis'  percentage  points  are  then 
used  as  a  standard  by  which  to  compare  approximations. 
The  omitted  percentage  points  for  approximation  3  are 
those  which  occurred  in  cases  v;here  V2  £  (ni  -  1)  v  i  +  1 . 

4 . 3      Summary 

A   comparison  of  the  preceding  tables  reveals  the 
following  generality  for  all  approximations.   The  S%   points 
are  uniformly  more  accurate  than  the  corresponding  1% 
points,  and  the  accuracy  of  both  the  5%  and  1%  points 
increased  with  Va  ,  while  decreasing  Vv'ith  both  vi  and  m. 

Seme  unexpected  results  are  also  revealed.   The 
51  and  1^  points  of  T  computed  from  the  x^  approximation 
(approximation  1)  are  surprisingly  close  to  the  exact 
X'alues  for  m  =  2  and  the  values  given  by  Davis  in  the 
cases  m  =  3  and  4,  even  for  small  vz-       In  comparing  this 
approximation  with  Pillai's,  the  y^  approximation  is 
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generally  more  accurate  for  both  the  5"6  and  1-6   points. 
Some  exceptions  occur  in  the  1%  points  when  both  vi  and 
V2  are  large. 

Although  it  was  hoped  that  approximation  2  would 
be  relatively  accurate  overall,  the  results  w'ere  rather 
disappointing.   This  approximation  is,  without  exception, 
less  accurate  than  the  others.   It  is,  however,  very 
simple,  and  fairly  accurate  for  large  values  of  V2. 

The  results  of  approximation  3  were  very  satisfying. 
Both  the  5"o  and  Iv    points  for  this  approximation  were  gener- 
ally more  accurate  than  Pillai's,  and  markedly  so  for  small 
V2.   The  decrease  in  accuracy  for  increasing  Vi  seems  to 
be  less  for  approximation  3  than  for  the  other  approximations 
The  fact  that  this  approximation  cannot  be  applied  when 
V2  £  (m  -  1)  V 1  +  1  should  not  be  too  disturbing,  since  in 
most  practical  situations  we  would  expect  V2  to  be  large, 
relative  to  v i . 

Although  the  comparison  of  percentage  points  is 
only  for  the  cases  m  =  2,  3  and  4  and  selected  values  of 
Vi  and  V2,  it  is  suspected  that  the  relative  accuracy  of 
the  approximations  would  be  maintained  for  general  m,  Vj 
and  V2.  It  is  therefore  suggested  that  the  distribution 
of  T  be  approximated  by 


T  *  ^'^'^         F 

.  V2  -m+i  '■v,V2-m+l 


' 

m 

V 

.1   h 

mv  1 

1 

1  =  1 

(4.3.1) 
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when  V2  >  (ra 
given  by 


1)  Vi  +  1.   In  this  case  the  value  o£  v  is 


V  =   nivi  (v2-m) 

V2-"()T1-1}  Vi+1 


(4.3.2) 


If  the  situation  should  arise  when  V2  £  (m  -  1)  vi  +  1, 
the  x^  approximation  is  suggested. 


in  which 


T  *  ax. 


in 


(4.3.5) 


a  = 


(V2 -1)  (vi+V2-ni-l) 
(v2-in)  (v2-ni-l)  (V2-in-3) 


(4.3.4) 


and 


V  = 


mv  1  (v2~Jb)  (v2-m-5) 
(V2  -1)  (V1+V2 -m-1) 


(4.3.5) 


CHAPTER  V 
THE    POWER   OF   HOTELLING'S    GENEK-XLIZED   T^^ 

5.1      Intro  duo  ti  071 

In  Chapter  I  it  was  noted  that  the  power  function 
o£  T  ^,  or  equivalently  of  T, 

6CA;in,vi  ,V2:t)  =  P(T  >  tlA),  (5.1.1) 

where  A  is  the  m  x  m  diagonal  matrix  of  elements  Ai,**",X^ 

defined  in  (1.1.23),  is  a  function  only  of  m,  Vj,  Va  and 

Xi  •••, X  .   This  power  function  is  known  to  be  an  increasing 
'    '  m        ^ 

function  of  each  \.     [5],  and  T^^  is  known  to  be  an  admissi- 
ble test  of  the  general  linear  hypothesis  [10]  .   Aside  from 
these  properties,  little  is  known  concerning  the  behavior 
of  the  power  function  of  T^  in  tlie  general  case. 

Pillai  and  Jayachandran  [30]  have  compared  the 
powers  of  four  test  criteria  (see  (1.1.24),  (1.1.25), 
(1.1. -26)  and  (1.1.27)]  for  testing  the  general  linear 
hypothesis  in  the  case  m  =  2  for  selected  values  of  vi 
and  V2.   On  the  basis  of  these  po\-;er  comparisons,  Pillai 
and  Jayachandran  noted  some  apparent  properties  of  these 
criteria,  summarized  below. 

1,   Roy's  largest  root  criterion  seems  to  have  the 

smallest  power  overall  of  the  four  test  criteria 
considered. 


92  - 


-  93  - 


2.  The  power  of  Pillai's  criteria  is  generally 
somewhat  greater  tJian  that  of  the  other  three 
criteria.   Some  exceptions  were  noted,  however, 
with  respect  to  the  likelihood  ratio  test  and 
Tq2  for  large  V2  when  Xi  and  X2  are  not  nearly 
equal . 

3.  Tq^  seems  to  have  greater  power  than  the  likeli- 
hood ratio  for  large  V2  except  when  Ai  and  X2 
are  nearly  equal. 

4.  None  of  the  four  criteria  considered  seem  to  have 
a  monotonic  pov;er  function  in  the  sum  Xi  +  X?. 
For  fixed  Xi  +  X2,  as  Xi  and  X2  tend  to  be  equal, 
the  powers  of  Pillai's  criteria  and  the  likeli- 
hood ratio  test  increase,  while  the  powers  of 
Tq2  and  Roy's  largest  root  decrease. 

With  reference  to  remark  4  above,  it  should  be 

pointed  out  that  for  any  two  sets  of  alternatives,  Ai  and 

A'z  ,    where  A.  :   i  =  1 , 2  ,  is  an  m  x  m  diagonal  matrix  with 

diafTonal  elements  X..:   j  =  l,'**,m,  such  that 


tr  Ai  =  tr  A2 , 


(5.1.2) 


T  ^  cannot  have  greater  power  against  one    of  these  alter- 
natives for  all  levels  of  significance.   In  order  to  see 
this ,  consider 


P(T  >  tJAi)  dt  =  ff(T|Ai) 


mv 1  +  t  r  A : 
V2  -  m  -  T 


(5.1.3) 


Equations  (5.1.2)  and  (5.1,3)  imply  that 


P(T  >  tiAi)  dt 


P(T  >  tJAz)  dt, 


(5.1.4) 
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so  that  it  cannot  be  true  that 

P(T  >  t|A^)  >  P(T  >  tlA.)  ,  (5.1.5) 

for  all  t  ^  0. 

Rather  than  give  further  power  coir.parisons  of 
various  test  criteria,  we  present  here  soine  integral 
bounds  on  the  power  function  of  T^^.   These  are  based  on 
a  result  due  to  IVald  [35]  and  Hsu  [16],  and  the  proper- 
ties of  completely  monotone  functions  [9].   In  addition, 
we  offer  an  alternative  proof  of  the  result  obtained  by 
Das  Gupta,  Mudholkar  and  Anderson. 

5.2      An    Upper   Bound   on    the    Integral    of   the 
Power   Function   of   Tq2 

In  1942,  l\'ald   [35]  and  in  1945,  Hsu  [16]  proved 
the  following  result  in  connection  with  the  analysis  of 
variance.   Suppose  we  are  given  zi,'*«,z   and  r^  with  the 
variates  independently  distributed  such  that 

2.  ^  N(c.,a2]:    i  =  1,---,P  (5.2.1) 

and 

r^  -  c^X  '(0).  (5.2.2) 

It  is  desired  to  test  the  hypothesis  that 

P  ,    =    ...    =    I      =    0  (5.2.3) 
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against  the  general  alternative.   Let  co*  be  the  critical 
region  o£  the  form 


w' 


(zi  ,• • • .z^,r2) 


i  =  l 


;.Vr'  >  c 
1  ' 


(5.2.4) 


where  c  is  chosen  so  that  oj-  is  of  some  desired  size,  a. 
The  power  function  of  u*, 


0 (n ;p,n)  =  P 


I   z  Vr^  >  cln 
i  =  l   ^ 


(5.2.5) 


depends  only  upon  p,  n  and  the  parameter 


1   ? 
h   .1    h' 


(5.2.6) 


Suppose  oj  is  any  other  critical  region  of  the  same  size  as 
u*  ,  and,  with 


i  =  (Ci,---,Cp)', 


(5.2.7) 


denote    its   r-ower    function   bv 


P^Ci;p,n) 


(5.2.8) 


For  fixed  t]  ,    let  S(_K,a,]])    denote  the  sphere  defined  by 


(5.  2. 6  J  so  that 


Arc,a,n)  = 


dK 


(5.2.9) 


Sa,o,ri) 


is    the    surface    area   of   tliis    sphere.      Then 
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2°^'T'P'^^^  1  AU,g,nJ  . 


S(C,a,n) 


6^(C;p,n)  d|_ 


(5.2.10) 


V/e  use  this  result  in  order  to  obtain  an  upper 
bound  on  the  integral  of  the  pouter  function  of  T^^ .   We 
have  an  m  x  vi  :natrix,  Y,  and  an  m  x  m  matrix,  S^,  ,  where 
Y  and  S2  are  independently  distributed  with 


and 


Y  -  N     rH,I  0  I) 


S2  -  r./^(l,v2,(0)] 


(5.2.11) 


(5.2.12) 


(see  (3..1.6)  and  (3.1.10)].   It  is  desired  to  test  the 
hvt:othesis 


against 


H  -  =(0) 
o 


H^:H^(0). 


(5.2.13) 


With   S2    distributed   as    in    (5.2.12)    we   may  write 
(see    [1]    or    [32]) 


S2    =    RR', 


(5.2.14) 


where  R  is  the  unique  m  x  m  lower  triangular  matrix  defined 
by  (5.2.14)  such  that 


r. 


and 


.  .  -  Y    .^,  (0)  :    i  =  I,-"  ,m 
11    '^V2-i  +  l 

N(0,1)  :        i  >  j 


r.  - 


(5.2.15) 


(5.2.16) 
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with  all  variates  independent.   Then 


tr  S2  =  tr  RR 


m   1 


I    I  r 
i=l  j=l  ^^ 


^mv2 


(5.2.17) 


From  (5.2.11) 


y.j  -  N(e..,l):    i  =    l,---,m;  j  =  l,.-.,vi.   (5.2.18) 

Let  CO*  be  the  critical  region  o£  size  a  for  testing 
(5.2.13)  of  the  form 


m   V 1 


S2  >  c 


(Y,S2):  I        I      Yi.Vtr 
i  =  l  j=l    ■' 

The  power  function  of  w*, 

So (n  ;m,Vi ,V2 :c) , 

depends  only  upon  m,  Vi,  V2  and  the  parameter 

m   V 1 


i=l  3=1   ^J 


(5.2.19) 


(5.2.20) 


=  tr  E; 


=  tr  A 


(5.2.21) 


[see  equations  (3.1.7)  and  (3.3.4)).   Then  from  the  result 
due  to  Wald  and  Hsu  (5.2.10) 


9S 


Bo(n;m,vi,V2:c)  1  AiA;l,v,,-)    \ 


S(A;n,Vi ,H) 
X  3(A;ra,vi ,V2:t)  dE ,  (5.2,22) 

where  S  (A;in,v  i  ,5)  denotes  the  sphere  defined  by  (5.2.21) 
and 


A(A;ni,Vi,H)  =  ]  dH ,  (5.2.23) 

JS(A;m,vi,H) 

its  surface  area.   From  (5.2.11),  (5.2.17)  and  (5.2.19) 


3o(n;in,vi,V2:c)  =  ^[x^^.Htr   -^)/x^^,HO)    >  c)  .   (5.2.24) 
Thus 

Pfx    ^(tr  A)/x    HO)  >  c] 

^  \(\-m^v,    -I  I  6(A;m,Vi,V2:t)  dH   (5.2.25) 

is  an  upper  bound  on  the  integral  of  the  power  function  of 
Tq^.   It  is  noted,  of  course,  that  c  and  t  are  chosen  to 
satisfy 

pfx   HO)/y        HO)    >  c] 

=  3[(0)  ;:n,Vi,V2:t] 

=  a.  (5.2.26) 

for  sone  level  of  significance,  a. 
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5.3      A    Lower  Bound,  on    the    Integral    of    the 

Power   Function    of   T  ^ 

■'       o 

The  lower  bound  on  the  integral  of  3  (A ;m,v i , V2 : t) 
is  based  on  the  properties  o£  completely  monotone  functions 
as  described  by  Feller  [9].   For  reference,  we  note  some 
of  the  basic  properties  of  these  functions. 

A  function,  $f6),  0  <  9  <  ~,  is  said  to  be 
completely  monotone  (6)  if  0(9)  possesses  derivatives  of 
all  orders,  ^    (9) ,  and 

(-1)^  $f")(e)  >  0;    0  <  G  <  -.  (5.3.1) 

Some  properties  of  completely  monotone  functions,  which 
follow  directly  from  the  above  definition,  are  listed 
below. 

1.  If  $(9)  is  completely  monotone  (9),  then 
(-1)  ^^       (0)  is  completely  m.onotone  (6). 

2.  If  <p(6)  and  T(9)  are  completely  monotone  (0), 
then 

Ri  (9)  =  $(9)  +  ^^(e)  (5.5.2) 

and 

R2(9)  =  $(e)  H-O)  (5.3.3) 

are  completely  monotone  (9). 

3.  If  5'(9)  is  completely  monotone  (8)  and  '•F(9)  is  a 
positive  function  v/hose  first  derivative  is 
completely  monotone  (9) ,  then 
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RsCQ)  =  ^?[^(9)]  (5.3.4) 

is  completely  monotone  (0) . 
The  main  property  o£  completely  monotone  functions  is  the 
follov/ing  theorem,  proved  in  Feller  [9]. 

Theorem:      A  function  $(6)  on  0  <  9  <  <»  is  completely 
monotone  (9)  if  and  only  if 


4(9) 


f  00 


J  0 


e'^^  FCdx) 


(5.3.5) 


where  F(x)  is  a  (possibly  infinite)  measure  on  [0,°°] 
It  follows  that  if  *(9)  can  be  written  in  the  form 


$(0) 


h(x)  e'^^  dx, 


(5.3.6) 


then  '±'(0)  is  completely  monotone  (9)  if  and  only  if  h(x) 
is  nonnegative  on    [0,°°]. 

Recall  the  representation  of  T  given  in  (3.2.34) 


T  =  I 


m      u.  ■ 

1 

1  =  1   1 


(5.3.7) 


where,  conditional  on  the  m  x  m  random  orthogonal  matrix, 
K  (see  equations  (3.2.25),  (3.2.29),  (3.3.30)  and  (3.3.31)) 
ui^,***,u  ^  are  independently  distributed  with 


and 


u-MK  -  X^/  (y^]  :    i  =  I,---  ,m 


y-  =  (KAK')ii 
m 


(5.3.8) 


(5.3.9) 
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The  variates  di,  •••,(!   are  the  latent  roots  of  S2  .   For 


convenience,  write 


a. 


i  ■■=    1 ,  •  •  •  ,m, 


(5,3.10) 


Thus,  using  (.5.3.7),  (5.3.8),  (5.3.9)  and  (5.3.10),  for 
fixed  a'  =  (ai,*'',a  )  and  K,  we  have 


m 


T|a,K  ^  I      ax,,  '  (y.) 
i=l   ^  ^1    ^ 

We  wish  to  consider  the  function 


(5.3.11) 


$(6)  = 


P(T  >  t|a,K,A)  dt 


e'^^  ds.    (5.3.12) 


The  inside  integral  of  (5.3.12)  is  nonnegative  and  hence 
from  (5.3.6),  'I'(G)  is  completely  monotone  (9).   Integration 
of  (5.3.12)  along  with  (5.3.11)  yields 


m 


(6)  =.  1  I      a.(v,  -  Yi)  -  ^ 

1=1  ° 


m 


+  —  n   (l  +  2ea.)   •'   exp 


-Vi/2 
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i  =  l 


ea.Yi 


1  +  2ect: 


Writing 


g-(6)  =  (1  +  25)-^i/2  g^p 


OYi 
1  +  2Z 


(5.3.13) 


(5.3.14) 


we  have 
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HQ)    =   ^     I      c..(v,  -  y.)  -  ^.  gL  _n   g.(9a.). 

(5.3.15) 

It  should  be  noted  that  gj_C5)  defined  in  (5.3.14)  is  the 
Laplace  transform  of  a  randon  variable  distributed  as 
X   ^(y-).   Hence,  fron  (5.3.6),  g^(5)  is  coiripletely 

monotone  (o)  . 

A  function,  f  (xi  ,  •  •  •  ,x^)  ,  is  convex  in  Xi,---,x^^ 
if  and  only  if  its  Hessian,  the  m  x  n  matrix,  H,  of 
second-order  partial  derivatives, 

h..  =  ^-^^i — :   i  =  l,---,Tn;   i  =  l,---,n,     (5.3.16) 

Ij     oX^cX^. 

is  nonnegative  definite.   Consider  the  second-order  deriva- 
tives of  5(9)  v;ith  respect  to  ^i,  *"•.%•   Direct  differ- 
entiation shov."s  that 


i!:ne)    ^^^  \"^^^;    ^     'i  ^    "^   n   g-Oa.)    (5.3.17) 
3a. 2  (1  +  ^a^9)2 j=l   -    " 


and 


illClI  =  J!Li:^ii^!ilV    ^   grg.).     (5.3.18) 
3a. 3a.     1  +  29a.  J  (1  +  Z6a.)  ,_-,  =k  ^   k^ 
13  1  3   ^--l- 


Thus  the  Hessian  becomes 

m 


H  =  Hi   n   g. (9a.) ,  (5.3.19) 

i  =  l   ^ 
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where  Hi  is  the  tii  x  m  matrix 


Hi  = 


(vi+Yi)  (vi+Yj^^) 
(l+29aiJ (l+29a  J 


»    •    • 


(l  +  26aiJ(i  +  2eaJ- 


(vi+Y-^)(vi+Yj„  +  2) 
(l42ea  }2 


The  matrix  Hi  can  be  written  as  the  sum 


(5.3.20) 


Hi  =  H2  +  H3, 


(5.3.21) 


where 


H 


2  - 


Cv  1  +Y  1  ) 


ri+2Bai J 


(vi-Yj 


(1  +  290.  J 
^     m 


(vi+Yi) 
(i+29ai) 


and 


i^^^yJ  1 


m- 


'  (1+Z9a  ) 
^     m 


(5.3.22) 


H3  = 


2(vi+Yi) 
(l+2Saij2 


2(vi+Y^^ 

(l  +  29ci  j2 
^     m 


(5.3.23) 
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Since  H2  is  nonnegative  definite  and  H3  is  positive  definite, 
Hi,  and  therefore  the  Hessian  H,  is  positive  definite. 
Using  the  definition  of  convexity,  $(9)  given  by  (5.3.12) 


or  (5.3.15),  is  convex  in  ai,-**. 


a 


m 


Jensen's  inequality  states  that 
E(f(xi,---,x^))  >  f(F(xa),---,ff(xj)       (5.3.24) 

whenever  f(xi,---,xj  is  convex  in  Xi,--,x^,  provided  of 
course  that  these  expectations  exist.   From  (5.3.10)  and 
(3.3.18) 


^(«i^  =  ^T^m^'-         ^    ^    l'"-'^- 


(5.3.25) 


Replacing  a-  by  its  expectation  in  (5.3.15),  denoting  this 
new  function  by  $*(e),  and  applying  Jensen's  inequality  to 
$(9) ,  we  have 


where 


£($(9)]  -  ^^(9)  >  0, 
a 


.  0 


P 


E(T)  >  t|K,/V 
a 


dt 


-  6s  , 
e     ds 


(5.3.26) 


m 


m 


i  =  l  /  -  2.  .  J-  n   g. 


V2 


-^--J 


(5.3.27) 
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Letting 


-FCe)  =  ^(^9))  -  5^(9), 


(5.3.28) 


we  wish  to  show  that  'F(G)  is  completely  monotone  (9).   We 
need  to  show  that 

(-1)">F^^^(9)  >  0:   n  =  1,2,---.         C5.3.29) 


Now 


(-1)^4''^"^  (9)  =  (-1)^  ^ 


d9 


g[^(9))  -  <5aC9) 


I 


=  E 
a 


(-1)^$^^^(9) 


(-1)"$^'^'^^  (9)  .  (5.3.30) 


Since  (-1)"$^  ^^^  (9)  is  simply  (-1)^^$  ^^^  (e)  with  ai,---, 
replaced  by  their  expectations,  (5.3.29)  will  be  proved 
if  we  can  show  that  (-I)"'? '^^'^  (9)  is  convex  in  ai,''',a 


a 


m 


and  applying  Jensen's  inequality.   Consider 


3a^2 


(-1)"$'^^^(9) 


(-1) 


iH 


n  d'^   3^4(9) 


j^n   3a. ^ 
d9      1 


(5.3.31) 


and 


9a. 3a . 
1  3 


(-l)^c5^^^(9) 


1 


=  (-1) 


n 


in 


d'^  dH(_e) 


inH  3a.  3a. 
d9     X         J 


(5.3.32) 


From  (5.3.19),  (5.3.31)  and  (5.3.32),  the  Hessian  of 


,n 


'H*  =  (-1)^  L^   H 


de 


n 


>v 
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.       d9 


m 


n 


Hi      n      g.(9a.] 
i  =  l      ^        ^ 


V 


n     d  rr 


(-1)-   ^-  Hi 


de 


n 


+    Hi 


.-1    d 


n        m 


(-1)"   ~n      n    g.(8a.) 

de^  i=i  ^ 


(5.3.33) 

where   Hj    is    the   positiAre    definite   matrix   of    (5.3.20).      From 
(5.3.21),    (5.3.22)    and    (5.3.23), 


(-1)"   ^H,    =    (-1)^   ^.    (H,    .   H3] 

de  de 


=  H4  +  H4 '  +  H5, 


(5.3.34) 


v\?here 


H 


k    - 


(vi+Yi)n!  (2ai) 


n 


(l  +  29ai)- 


(vi+Yj,pn!  (2a^) 


n 


(l-^29a^) 
^  m 


n+1 


(v.^Y.)  ...  ^^^y' 

(l  +  29ai)'  '(l  +  20aj^) 


(5.3.35) 


and 
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H. 


2(vi+Yi)(n+])!  (2ai)" 


Cl+2eai) 


n+2 


2(vi+Yj„)Cn+l)!(2a^)" 


(l+2eai) 


nTT 


(5.3.36) 

so  that  C-1)"  d^/de'^  Hj  is  positive  definite.   Recall  that 
g.  (6)  is  defined  in  (5.3.14)  is  completely  monotone  (6). 
The  function 


h^(9)  =  ea^ 


(5.3.37) 


is  a  positive  function  whose  first  derivative  is  completely 
monotone  (9) .   Using  property  3  of  completely  monotone 
functions  (5.3.4),  we  see  that  g..  (9a.)  is  completely  mono- 
tone (9).   Repeated  application  of  property  2  (5.3.3) 
implies  that 


m 


G(9)  =   n   g-(ea.) 
i  =  l   ^    ^ 

is  completely  monotone  (9) .   Thus 


(5.3.38) 


jH 


de 


(5.3.39) 
and  hence,  from  (5.3.33)  the  Hessian,  H* ,  of  (-1)  <f>^'-^(e) 
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is  nonneqative  definite.   This  shows  that  (-1)^$ '-^-' (9)  i: 


convex  in  ai,'**,aj^,  v/hich  in  turn  shows  that 


"/(e)  =  g($(e)]  -  $^(6) 


(5.3.40) 


is  completely  monotone  [9) . 

The  definition  of  <5(9)  is  given  in  (5.3.12)  and 
our  definition  of  *  (9)  given  in  (5.3.27)  yield 


^(Q)    = 


f(P(T  >  t|a,K,A))  dt 
is  £ 


P(^(T)  >  t|K,A)  dt 


a 
s   — 


e     di . 


(5.3.41) 


Equation  (5.3.6)  and  the  complete  monotonicity  of  "Hid) 
imply  that 


s  ^ 


£'fP(T  >  t|a,K,A)]  dt 


J  s 


P(g(T)  >  tlK,A]  dt, 


(5.3.42) 


for  all  s,  0  £  s  <  ^^.   From  (5.3.7),  (5.3.8),  (5.3.9) 
(5.3.10)  and  (5.3.25)  and  since  K  is  orthogonal. 


p(g(T)  >  t!K,A)=  P 


m 


y     a .  u .  ' 
i  =  l   ^  ' 


>    t!K,A 
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V2 


^^—^  y         Htr  A)  >  t 
-m-1  '^'mvi  ^     ^ 


=    P 


Xj^^^ntr  A)  >  (V2  -  m  -  1)  t 


(5.3.43) 


The  last  equation  follows  from  the  fact  that  the  one 
immediately  preceding  this  does  net  depend  on  K.   Equations 
(5.3.42)  and  (5.3.43)  then  imply  that 


r°o 


^(P(T  >  t|a,K,A)}  dt 
s  - 


,00 


P(X^^/  (tr  A)  >  (V2  -  m  -  1)  t)  dt , 


(5.3.44) 


for  all  s,  0  _<  s  <  =».  Since  (5.3.44)  is  true  for  all 
orthogonal  K,  it  is  true  unconditionally  on  K.  This, 
along  with  the  fact  that 


g(P(T  >  t|a,A))  -  P(T  >  tiA), 


(5.5.45) 


yields  a  lower  bound  on  the  integral  of  tlie  power  function 


of  T  % 
o  ' 


S  (A;m,\)i  ,V2  :  t)  dt 


-  r  ^^V./^'^^  A)  >  (V2  -  m  -  l)t)  dt, 
■'  s 

(5.3.46) 


for  all  s ,  0  <  s  <  «>. 
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5.4      The   Monotonicitv    of   fm    Power 


Function    of  T    ^ 


Let  A : (m  x  m)  and  A* : (m  x  m)  be  diagonal  matrices 
with  diac^onal  elements  Xi,'**,X„  and  Xi*  ,•••  ,\„=^  respec- 
tively,  and  suppose  that 


X .    >    \ .*      for  some  i , 
1—1  ' 

and  (5.4.1) 

\.    =    \.*      for  i  ^  i. 
3  3  ^ 

We  shall  demonstrate  that 

BCA;m,vi  ,V2  :t)  ^  B  (A*  ;m  ,v  i  ,\;2  :  t)  (5.4.2) 

for  all  t  >_   0. 

With  the  representation  of  T  given  in  (3.2.34)  and 
an  obvious  notational  convention,  we  have 

m     u.  ^ 
T    =      I      ^  (5.4.5) 

i  =  l   i 


and 


u.  ''^ 

T=*  =  I    -^,  (5.4.4) 

i 


.n  which    [see    equations    (3.2.28)    and    (3.3.6)) 


and 


u.2|K    -    X^/(Yi)  (5.4.5) 


Ui*MK    -    X^    '(Yi*),       ■  (5.4.6) 


-  Ill 


where    [see   equation    (3.3.5)) 


m 


Y 


.  =   y   k.  .2  A. 

1      iii    ^J      3 


and 


C5.4.7) 


m 


Y. 


A  - 


y   k.  .2  X.*. 
j=i    ^J     ^ 


(5.4.8) 


Recall  that  conditional  on  K,  Ui^,''«,u  ^  are  mutually 


m 


independent  as  are  Ui*^ ,•••, u  *^ .      With 

m 


^i  "  dT  •   i  =  1,- •• ,m, 
i 


(5.4.9) 


and  for  fixed  o.'    =    (ai,**«,a  ),  consider  the  function 
—  m 

Q(e)  =  [   [P(T  >  t[a,K,A)  -  P(T=^  >  t|a,K,A*))e'^'^  dt 


m 


I      a.u. 2  >  tla,K,A 
i  =  l   1  1 


P(I  »-i"i*'  >  t|a,K,A*) 


e     dt , 


(5.4.10) 


Equation  (5.4.2)  is  proven  if  we  can  shox'/  that  Q(9)  is 
completely  monotone  (9)  for  all  a-  >  0:   i  ===  l,''',m. 
Direct  integration  of  (5.4.10)  yields 

1 


Q(9)  = 


exp 


m 


9   y   a.u.*Ma,K 

-  ^T      11      '  — ' 
1  =  1 
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exp 


m 


e   y   a  .  u  .  M  a  ,  K 

1  =  1 


(5.4.11) 


Using  (5.4.5)  and  (5.4.6)  and  taking  expectations  condi- 
tional on  K,  we  have 


Q(0)  =  Q  R(e)  E 
K 


exp 


exp 


m 


m 


.1    h^^^  '<i 
1=1 


I  S.(e)y- 
i  =  l   ^    ^ 


a,K 


(5.4.12) 


where 


m 


and 


R(e)  =  n  (1  +  2ea.) 

i  =  l 


ea. 


•vi/2 


(5.4.13) 


^i"^^^  "  1  +  Ida. 
For   r  =  1 , • • • ,m   let 


(5.4.14) 


q  (6)  =  F 


exp 


r  m 

I      6.(9)  Y.*  -   I    5.(6)  y. 
i=l   ^     ^    i=r+l   ^     ^ 


exp 


r-1  m 

I      6.(6)  Y-*  -  J   S.(e)Yi 
1=1  i=r 


a,K 


(5.4.15) 


Then  v.'e  mav  write 
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m 

V 


(5.4.16) 


Q(e)    =   ^  R(9)    I      q    (0). 
r=l 

From  equations    (5.4.1),    (5.4.7)    and    (5.4.8), 


>    P[y/    >    slYi*,-..,Y^.;L*'^r  +  l'*"'^m^       (5.4.17) 

for  all  s  >  0  and  r  =  l,««',Tn.   Therefore  [see  equation 
(5.3.6)] 


^^(5)  = 


P(y  >  s\':,*,"',y^.i*,y^^i,"',yj 


-  P(Yj.*  >  s|y 


1  '    '  ^r-1  '  ^r  +  1 '    '  'm-' 


-5s  . 
e     ds 


<S  K 


-6y/ 
e    ~ 


-5Y. 


I. 


v"*»»v        V       •••V 

11  >    ''r-1    r+1'     iTi 

(5.4.18) 


is  completely  monotone  (6).  Since  3- (9)  (see  (5.4.14))  is 
a  positive  function  with  a  completely  monotone  derivative, 
then  from  (5.3.4)  , 


1    J   -6r'^Q)Yr* 


,1;  [3  (9)1  =  TT-T-TTT-  i'   e 
^r^  1-^  J  J    B^TT  K 


-Br(B)Yr 


Yi  ,    ,f^.;^  ''r  +  1'    '  ^m 


(5.4.19) 
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is  completely  monotone  (e) .   This  in  turn  implies  that 
r  r 


I  i^A'»  -  rAmr-  t^-^-") 


is  completely  monotone  (6) .   Now 

a 

T 

r 

is  completely  monotone  (9),  and  hence,  using  (5.3.3)  and 
(5.4.20),  ^q,(e)  is  completely  monotone  (9).   Direct 
differentiation  shows  that  R(e)  (see  equation  (5.4.13)) 
is  completely  monotone  (9)  and,  using  (5.4.16)  along  with 
an  application  of  property  2  of  completely  monotone 
functions  (see  equations  (5.3.2)  and  (5.3.3)),  implies 
that  Q(9)  is  completely  monotone  (9).   Equation  (5.3.6) 
implies  that 

P(T  >  t|a,A)  >  P(T*  >  tIo.,A*)  (5.4.22) 

for  all  t  >  0.   Since  (5.4.22)  is  true  for  all  a,  it  is 
true  unconditionally.   This  demonstrates  the  result  of 
Das  Gupta,  et   al .  , 

6(A;m,vi  ,V2  :t)  >  B  (A*  ;in,  vi  ,V2 :  t)         (5.4.23) 
for  all  t  ^  0. 

5. 5  Summary 

Little  is  known  concerning  the  behaAdor  of  the 
power  function  of  T  ^  for  general  m,  vi  and  V2,  with  the 
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exception  that  it  is  an  increasing  function  of  each  X 
Upper  and  lo'.ver  bounds  on  the  integral  of  the  power 
function  of  T  '^  vere    derived  in  Sections  5.2  and  5.3. 
These  are 


S(A;m,,vi  ,  5) 


6 (A;m,vi ,V2 : t)  ds 


<  ACA;ni,Vi  ,  5)  P 


X„  ,  ^(tr  A)/x  ^   2(0)  >  c 


(5.5.1) 


and 


3(A;m,Vi ,V2 :t)  dt 


PfXj^^^2(tr  A)  >  (V2  -  m  -  l)t 


dt 


(5.5.2) 


where  in  (5.5.1)  t  and  c  are  chosen  such  that 


3[(0) ;m,vi ,V2:t)  =  P 


X    '(0)/x^   '(0)  >  c 


=  a 


(5.5.3) 


for  some  desired  level  of  significance,  a. 

Although  we  have  not  yet  been  able  to  prove  this, 
it  xvould  seem  most  likely  that  g  (A  ;m,  v  i  ,  V2  :  t)  is  an  in- 
creasing function  of  vz-      If  this  is  true,  an  upper  bound 


can  be  placed  on  the  power  of  T   ,  since 
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e(A;m,Vi,V2:t)  =  P(T^'  >  tJA)  ' 

=  PCv2  trY'Sz'^Y   >  t|A) 

=  P(tr  Y'CSo/vz)-^  Y   >  t|A]  . 

(5.5.4) 

If  the  conjecture  is  true,  then 

8CA;m,Vi,V2:t)  <  lin   Pftr  Y'CSa/va)"'  Y   >  1 1  A) 

=    P[x^^^^Htr   A)    >    t'),  (5.5.5) 

i\[here    t'    sa.ti3fies 

S[(0);in,Vi,V2:t)    =    P(x^^/CO)    >    t'] 

=   a.  (5.5.6) 

We  reneat  that  this  is  merely  a  conjecture,  but  would  seem 
a  fruitful  area  for  future  research. 
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